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description  of  available  atmospheric  models.  Trajectory 
prediction  programs  are  reviewed  and  a  revised  program  is  used 
to  conduct  a  simulation  of  balloon  trajectories.  Balloon 
locations  at  fixed  times  are  analyzed  for  variability.  The 
study  concludes  that  high  altitude  balloons  have  some 
potential  for  use  as  surveillance  platforms  for  limited 
periods  of  time. 
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I .  INTRODUCTION 


A.  PURPOSE 

The  Persian  Gulf  War  of  1990-1991  included  the  firing  of 
more  than  80  Scud  missiles  against  Coalition  forces  and 
Israel.  This  was  actually  the  second  employment  of  theater 
ballistic  missiles  against  United  States'  forces;  Germany's 
use  of  V-1  and  V-2  rockets  in  World  War  II  was  the  first. 
While  the  military  importance  of  mobile  tactical  ballistic 
missiles  can  be  argued,  we  learned  from  the  Persian  Gulf  War 
that  such  missiles  are  a  threat  and  they  are  neither  easily 
detected  nor  easily  destroyed.  The  purpose  of  this  work  is  to 
examine  using  high  altitude  balloon  systems  with  sensor 
payloads  to  increase  our  ability  to  detect  and  to  locate  enemy 
mobile  tactical  ballistic  missile  systems. 

B.  RESEARCH  QUESTIONS  AND  LIMITATIONS 

The  basic  hypothesis  of  this  study  is  that  employing 
balloons  as  surveillance  platforms  at  altitudes  in  excess  of 
21  kilometers  is  feasible.  Possible  research  questions 
include : 

•  How  will  balloons  drift  at  higher  altitudes? 

•  Will  either  "strategic"  global  circumnavigating  balloons 
or  "tactical"  theater-employed  balloons  pass  over  or 
nearly  over  a  designated  target  area  to  allow  some  set  of 
these  balloons  to  provide  continuous  coverage  of  a  target 


•  What  are  the  operational  requirements  of  a  balloon 
surveillance  system:  e.g.  how  many  systems  are  needed  for 
coverage  of  an  area  and  what  sensors  are  best? 

This  study  focuses  on  the  question  of  balloon  drift:  that 

is,  where  balloon  surveillance  systems  drifting  at  high 

altitude  would  go.  The  question  is  further  narrowed  to 

balloon  location  at  specific  times  after  launch:  ten  days,  one 

month,  and  one  year.  These  locations  were  used  for  a 

statistical  analysis  of  the  variance  in  mean  location  of 

drifting  high  altitude  balloons.  Thus,  the  perspective  is 

limited  by  these  discrete  and  disparate  time  points.  The 

research  presented  in  this  thesis  contains  no  classified 

information . 


ZI.  BACKGROUND:  PERSIAN  GOLF  WAR  SCUD  CAMPAIGN 

A.  SCUD  CAMPAIGN  EFFORTS 

When  Iraq  first  launched  its  Scud  missiles  during  the  war, 
conventional  wisdom  was  unanimous:  the  Scud  missiles  were 
militarily  insignificant.  The  reasons  were  that  Iraq  had  few 
missile  systems  and  those  they  had  were  inaccurate  and 
unreliable.  Unfortunately,  conventional  wisdom  concerning  the 
political  significance  of  the  scud  attacks  was  equally 
unanimous:  the  attacks  could  easily  bring  Israel  into  the  war, 
with  the  subsequent  result  that  the  Allied  Coalition  would 
disintegrate.  Given  this  possibility,  the  Allies  had  no 
choice  but  to  mount  a  major  effort  against  the  Scuds. 

It  is  reported  that  five  percent  of  the  air  sorties  flown 
during  the  war  were  against  the  Scuds  (Dunnigan,  1992,  p. 
155)  ,  but  that  figure  does  not  indicate  the  level  of  effort 
actually  directed  against  the  Scud  threat.  The  ratio  of  Scud 
missions  flown  out  of  the  total  number  of  strike  aircraft 
sorties  flown  would  be  much  more  than  five  percent.  Also, 
significant  non-air  assets  (Patriot  units,  launch  detection 
sensors  and  communication  assets)  were  employed  for  Scud 
defense.  With  this  level  of  effort.  Scud  operations  were 
degraded,  but  post-war  analysis  has  indicated  that  destroying 
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In  sum,  zero  pressure  balloons  can  be  used  to  reliably 


Scud  launchers  was  a  difficult  job  at  best  because  of  the 
complexities  involved  with  detection  and  engagement. 


B.  PERSIAN  GULF  SCUD  CAMPAIGN  DIFFICULTIES 

Initial  claims  reported  by  the  services  and  published  by 
the  news  media  during  and  immediately  after  the  war  suggested 
that  Allied  forces  were  extremely  effective  in  detecting  and 
destroying  mobile  launchers,  albeit  usually  after  launch. 
Upon  closer  examination,  however,  it  has  become  clear  that  we 
were  not  very  effective  in  doing  so.  In  fact,  "there  are  no 
confirmed  reports  of  any  Scud  launchers  that  were  destroyed  by 
Coalition  forces  during  the  war."  (Israel,  1993)  We 
particularly  had  difficulty  finding  the  launchers,  which  fired 
missiles,  quickly  enough  to  engage  them  before  they 
repositioned:  sensors  in  use  were  inadequate  to  provide  the 
range,  accuracy,  and  timeliness  needed.  Costs  and  operational 
characteristics  of  overhead  systems,  to  include  remotely 
piloted  vehicles  and  combat  aircraft,  precluded  continuous, 
low-cost,  low-risk  surveillance  of  an  area.  Shortly  before 
the  Persian  Gulf  War,  the  U.S.  Army  Strategic  Defense  Command 
defined  several  requirements  for  sensors  to  detect  and  to 
locate  mobile  missiles.  Future  sensors  must: 

•  Provide  broad  area  coverage. 

•  Provide  long  dwell-time  coverage, 

•  Provide  day/night/all-weather  coverage. 

•  Overfly  enemy  areas. 
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•  Identify  potential  targets  as  mobile  Scuds. 

•  Report  in  near-real  time. 

•  Cue  attack  assets.  (Severance,  1990,  pp.  28-29) 

Vfhat  systems  can  meet  these  requirements?  Our  Gulf  War 
experience  could  be  used  to  argue  that  the  U.S.  currently  has 
no  sensor  systems  that  satisfy  all  of  these  needs.  However, 
a  free-floating  high-altitude  balloon  surveillance  system  may 
be  able  to  do  so. 

C.  BALLOON  SURVEILLANCE  SCENARIO  1995 

Recent  Iraqi  behavior  on  several  fronts  has  raised  the 
possibility  of  military  confrontation.  Slaughtering  of 
Shiites  in  the  south,  heightened  tensions  directed  against  the 
Kurds  in  the  north,  Iraqi  defiance  of  United  Nations 
inspection  teams,  and  Iraqi  saber-rattling  about  gaining 
control  of  Kuwait,  long  considered  by  Iraq  as  it's  I9th 
province,  all  indicate  an  Iraqi  intention  to  break  out  from 
underneath  stifling  U.N.  sanctions  through  military  action. 
United  States  (U.S.)  forces  stationed  in  Saudi  Arabia  and 
Israel  are  prepared,  however,  with  limited  defensive  units  in 
place  and  with  alerted  strike  forces  on  board  ships  in  the 
Persian  Gulf.  Seeming  to  have  learned  from  his  failed  1990 
invasion  of  Kuwait,  Saddam  Hussein  daringly  launches  an  all- 
out  attack,  simultaneously  crossing  the  border  into  Kuwait 
with  ground  forces  and  launching  Scud  missiles  at  both  Saudi 
Arabia  and  Israel.  Air  raid  sirens  warn  friendly  personnel  to 
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take  cover  and  to  don  protective  masks.  Improved  Patriot 
anti-missile  rockets  streak  from  the  ground  to  intercept  the 
Scud  missiles  before  their  targets  can  be  reached.  The  dark 
night  sky  lights  up  brightly  from  the  fireworks  of  intercepts, 
and  the  Scud  missiles  are  obliterated  by  the  kinetic  energy 
impacts  high  in  the  atmosphere.  Friendly  personnel  breath  a 
little  easier  as  the  ''all-clear”  is  sounded  and  the 
counterattack  against  the  Scud  Transporter  Erector  Launchers 
(TELs)  begins.  Despite  an  inability  to  engage  the  TELs  before 
their  departure  from  the  launch  area,  U.S.  forces  are  able  to 
track  the  TELs  to  their  "hide"  locations  using  reconnaissance 
balloons.  Upon  the  launch  of  the  Scuds,  the  Infrared  sensor 
carried  in  the  balloon  payload  detects  a  launch  plume,  and 
then  triggers  a  Moving  Target  Indicator/Synthetic  Aperture 
Radar  (MTI/SAR)  to  point  to  the  launch  area.  The  MTI/SAR 
focuses  upon  the  vehicle  at  the  launch  point,  which  the  SAR 
processing  system  identifies  as  a  mobile  launcher.  The  MTI 
"locks"  onto  the  launcher,  tracking  it  as  it  repositions  from 
the  launch  location  to  a  hide  location.  That  location  is 
passed  to  an  in-theater  ground  station  which  directs  the 
mobilization  of  alert  aircraft  from  a  carrier  to  attack  the 
Scud.  Smart  munitions  are  subsequently  launched  to  penetrate 
a  sand-covered  bunker,  where  the  TEL  is  destroyed. 


III.  BALLOON  USE  AND  TECHNOLOGY 

A.  BALLOON  HISTORY 

1.  Buoyancy  and  Lift 

Archimedes  first  discovered  the  principle  of  buoyancy, 
that  a  body  immersed  in  fluid  is  buoyed  up  by  a  vertical  force 
equal  in  magnitude  to  the  weight  of  the  displaced  fluid,  in 
about  240  B.C.  The  object  of  a  balloon  is  to  displace  a  large 
weight  of  air,  thus  gaining  a  buoyant  force  equal  to  the 
weight  of  the  air  displaced.  Buoyancy  is  defined  by  the 
equation  Bouyancy=(p^^)  {Vol^)  (g)  (1-M,„/M^,)  =  (g)  (p^rP,J 

which  shows  that  balloon  lift  is  dependent  upon  the  volume  of 
the  balloon  (Vol,„)  and  the  difference  between  the  internal 
"lifting  gas"  density  (p,„)  and  the  external  air  density  (p^) , 
with  g  being  the  acceleration  due  to  gravity.  This  equation 
indicates  that  any  gas  having  a  molecular  weight  less  than  the 
molecular  weight  of  air  is  a  potentially  useful  lifting  gas. 
A  comparison  of  several  of  these  gases  and  their  lifting 
capability  is  shown  in  Table  1.  The  medium  with  the  smallest 
molecular  weight,  hydrogen,  provides  the  greatest  lift,  with 
helium  next  best.  (Morris,  1975,  pp.  IV-4  -  IV-7) 


* 


TABLE  1.  LIFT  GAS  COMPARISON 


Gas 

Molecular 

Lift  Index 

Weight* 

(1-M,/M.) 

Vacuum 

0.000 

1.000 

Hydrogen 

2.016 

0.930 

Helium 

4.003 

0.862 

Ammonia 

17.030 

0.412 

Water  Vapor 

18.010 

0.378 

(Steam) 

Air 

28.960 

0.000 

*  Standard 

Temperature  and 

Pressure  (STP) 

M,  =  Molecular  Weight  of  Gas 

M,  =  Molecular 

Weight  of  Air 

A  review  of  several  key  characteristics  of  a  few 
lifting  gases  reveals  why  helium  is  a  common  choice: 

•  a  vacuum  is  ideal  but  obviously  difficult  to  contain. 

•  hydrogen  is  readily  available  and  inexpensive,  but  it  is 
highly  flammable. 

•  helium  is  readily  available,  although  it  is  moderately 
expensive. 

•  ammonia  is  toxic  in  high  concentrations  and  it  liguifies 
easily. 

•  steam  liquifies  too  readily. 

•  air  is  inexpensive  but  it  provides  lift  only  if  the 
density  of  the  air  inside  the  balloon  is  less  than  the  air 
outside  the  balloon  (hot  air  balloons).  (Morris,  1975,  p. 
IV-5) 

Table  II  shows  examples  of  the  lifting  gas  required 
for  buoyancy  for  selected  payload  weights.  It  should  be  noted 
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TABLE  II.  HELIUM  NEEDED  FOR 
SELECTED  PAYLOAD  WEIGHTS 


Payload 

Balloon 

Lifting  Gas 

Weight 

Diameter 

Quantity 

(Ib/kg) 

(ft) 

(ft*) 

100/  45 

45 

1520 

500/225 

130 

7600 

1000/450 

200 

15200 

that  the  buoyancy  of  the  balloon  must  be  greater  than  the 
weight  of  the  payload  for  the  system  to  ascend:  the  balloon 
will  not  lift  a  payload  if  the  buoyant  force  is  less  than  the 
forces  holding  the  balloon  down.  Thus,  the  weight  of  the 
payload  suspended  below  the  balloon  determines  the  amount  of 
lifting  gas  (and  the  size  the  weight  of  the  balloon)  needed  to 
ensure  sufficient  buoyancy.  (Morris,  1975,  pp.  IV-A  -  IV-7) 
2.  Balloon  Invention  and  Development 

It  was  not  until  1250  A.D.  that  Roger  Bacon  suggested 
filling  a  vessel  with  some  fluid  lighter  than  air  to  apply 
Archimedes'  principle  in  the  atmosphere.  Cyrano  de  Bergerac 
proposed  using  a  belt  consisting  of  glass  phials  filled  with 
dew  to  help  him  ascend  to  the  stars  in  his  book  Voyages  to  the 
Moon  and  the  Sun,  published  between  1657  and  1662,  and  then 
Jacques  and  Joseph  Montgolfier  first  demonstrated  these 
concepts  in  November  1782.  Wondering  why  smoke  always  rose 
upwards  and  what  would  happen  if  the  hot  air  could  be 
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entrapped,  they  fashioned  a  lightweight  paper  bag  (their 
father  was  a  paper  manufacturer) ,  and  then  burned  paper 
underneath.  The  bag  rose  to  the  ceiling  of  their  room.  The 
Montgolfier's  repeated  their  success  outdoors,  and  then 
demonstrated  it  with  a  larger  balloon  before  the  Academy  of 
Sciences  in  Paris.  (Glines,  1965,  p.  5) 

Two  centuries  of  experimentation  in  shapes,  sizes, 
materials,  and  lifting  gases  have  followed.  Technological 
advancements  have  introduced  lighter,  thinner,  stronger,  and 
more  versatile  balloon  envelopes,  including  plastic  films. 
These  changes  have  resulted  in  increased  buoyancy,  decreased 
gas  permeability  and  longer  duration  flight.  As  a  result, 
today's  balloons  can  reliably  fly  higher  and  farther  than  ever 
before.  (Dollfus,  1961,  p.  97) 

3.  Military  Balloon  Use 

Balloons  were  first  used  for  military  purposes  during 
the  French  Revolution  in  1794  at  the  Battle  of  Fleurus,  where 
Frenchman  Jean  Coutelle  went  450  meters  aloft  in  a  tethered 
balloon  to  observe  enemy  formations  and  movements.  His 
observations  were  loudly  proclaimed  afterwards  as  key  to  the 
French  victory.  (Glines,  1965,  p.  98) 

Surveillance  balloons  were  successfully  used  during 
the  American  Civil  War  with  much  fanfare,  both  by  the  South 
and  the  North.  Union  use  was  much  more  successful,  led  by  the 
much-acclaimed  Professor  Thaddeus  Lowe  who  made  numerous 
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surveillance  flights  from  1861  to  1863.  Professor  Lowe's 
largest  contributions  were  made  during  the  Peninsular  Campaign 
in  support  of  General  Hooker.  Lowe  frequently  provided 
information  about  the  enemy  that  could  not  have  been  obtained 
without  his  Balloon  Corps  of  seven  balloons  and  a  Navy  vessel. 
Unfortunately,  neither  Professor  Lowe  nor  his  fellow  aeronauts 
were  ever  commissioned,  and  the  administration  of  his  Balloon 
Corps  was  passed  from  one  organization  to  another.  The  result 
was  that  the  Balloon  Corps'  transportation  assets  were 
frequently  taken  away  for  other  needs  so  it  became  impossible 
to  move  the  equipment:  Professor  Lowe  failed  to  reach  Antietam 
and  Gettysburg  quickly  enough  to  provide  observation  of  enemy 
movements.  The  difficulties  in  administrative  organization 
and  logistics  led  to  the  eventual  disintegration  of  the  Corps, 
with  Lowe  quitting  shortly  after  Gettysburg.  (Glines,  1965, 
pp.  101-109) 

The  Franco-Prussian  War  of  1870-71  saw  an  increase  in 
balloon  use  during  the  siege  of  Paris,  when  citizens  built  64 
balloons  which  they  used  to  escape  and  to  send  messages  over 
the  surrounding  German  troops.  By  the  end  of  the  war,  11  tons 
of  mail  and  164  people  had  been  successfully  airlifted  out. 
(Glines,  1965,  pp.  122-127) 


Balloons 

were  used 

during  both 

World 

Wars , 

predominantly  in 

the  form  of 

dirigibles,  or 

powered  air 

vehicles.  Count 

Ferdinand  von 

Zeppelin  had 

served 

with 

Professor  Lowe  during  the  American  Civil  War  and  he  had  become 


alarmed  by  the  developments  made  by  the  French,  who  had 
crafted  a  non-rigid  airship  which  could  fly  at  11  miles  per 
hour  powered  by  an  electric  engine.  The  non-rigid  French 
airship  could  provide  long-range  surveillance  and  carry  bombs. 
The  Count  built  a  new  rigid  airship,  named  after  Zeppelin 
himself.  The  Zeppelins  were  used  by  the  Germans  during  both 
World  Wars  for  resupply  and  bombing,  but  the  hydrogen-filled 
dirigibles  were  vulnerable  to  faster,  higher-flying  fighter 
aircraft  equipped  with  machine  guns,  so  their  use  was  limited. 
(Macksey,  1986,  pp.  48-49) 

Also  during  World  War  II,  the  Japanese  sent  9000 
parchment  balloons  carrying  incendiary  devices  toward  the 
United  States.  About  1000  balloons  reached  the  North  American 
continent,  but  the  only  recorded  casualties  were  five  children 
and  a  woman  at  a  picnic  in  Oregon.  The  U.S.  government 
requested  that  newspapers  not  report  that  the  balloons  were 
reaching  the  mainland,  and  the  media  blackout  led  the  Japanese 
to  conclude  that  their  program  was  not  worth  continuing, 
despite  the  fact  that  their  trans-ocean  success  rate  was 
actually  rather  good.  (Glines,  1965,  pp.  143-147) 

American  military  use  was  prompted  by  the  Japanese 
program's  success,  which  demonstrated  the  potential  of  using 
jetstream  winds  to  propel  balloons.  In  1950,  the  U.S.  Air 
Force  began  research  on  a  balloon  surveillance  system  to 
provide  reconnaissance  overflights  of  the  Soviet  Union.  This 
work  was  done  concurrently  with  the  development  of  the  U-2 


reconnaissance  aircraft  and  early  satellite  research.  In 
January  1956,  the  system  was  put  into  operation  with  balloons 
containing  cameras  launched  from  five  locations  in  Europe. 
They  drifted  across  Asia  at  13.6  kilometers  altitude,  taking 
pictures,  and  while  many  landed  over  foreign  territory,  about 
40  were  successfully  recovered  from  the  Pacific  Ocean.  The 
program  was  discontinued  on  March  1,  1956  following  Soviet 
protests  and  a  Washington  Post  story  on  February  10.  The 
program  was  able,  however,  to  photograph  over  1  million  sguare 
miles  of  the  Sino-Soviet  area,  at  a  cost  of  only  $48.49  per 
square  mile.  (Davies,  1988,  pp.  59-61) 

B.  CURREMT  BALLOON  TECHNOLOGY 

The  three  basic  types  of  balloons,  zero  pressure, 
superpressure,  and  sky  anchor,  are  shown  in  Figure  1  (after 
Lawrence  Livermore) .  Each  type  is  described  below. 
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Figure  l.  Basic  Balloon  Types 


1.  Z«ro  Prassur*  Balloons 

The  zero  pressure  balloon,  first  built  and  flown  by 
the  Montgolfier  brothers,  is  the  type  of  balloon  most  used 
today,  for  science,  military,  and  recreational  purposes.  The 
basic  principles  of  zero  pressure  balloons  have  not  changed  in 
the  200-plus  years  since  first  flight:  a  non-extensible  bag, 
open  at  the  bottom,  is  inflated  with  a  gas  (usually  hot  air, 
helium,  or  hydrogen)  lighter  than  the  surrounding  atmosphere. 
The  term  "zero  pressure"  is  used  because  the  internal  gas 
pressure  is  equal  to  the  external  gas  pressure  near  the  base 
of  the  balloon.  Slight  overpressure  is  maintained  higher 
inside  the  balloon  above  this  point  in  order  to  retain  the 
shape  of  the  envelope.  (Lawrence  Livermore,  1990,  p.  8) 

As  previously  noted,  balloon  buoyancy  results  from  the 
balloon's  volume,  or  atmospheric  displacement,  following 
Archimedes  principle.  From  the  equation  for  buoyancy  it  can 
be  seen  that  buoyancy  increases  as  the  lifting  gas  becomes 
less  dense;  and  that  buoyancy  decreases  as  the  lifting  gas 
becomes  more  dense  or  as  balloon  volume  decreases.  Since 
lifting  gas  pressure,  temperature,  and  density  are  related  by 
the  equation  P=(r^„)  (Pja,)  (®)  where  B  is  a  scaling  constant, 

the  altitude  of  a  zero  pressure  balloon  is  dependent  upon  the 
temperature  of  the  lifting  gas.  As  the  gas  temperature 
increases,  balloon  volume  increases  and  as  the  gas  temperature 
decreases,  balloon  volume  decreases.  In  the  first  case,  the 


balloon  becomes  more  buoyant  and  it  will  rise  in  altitude.  In 
the  second  case,  the  balloon  becomes  less  buoyant  and  it  will 
fall  in  altitude.  An  example  helps  illustrate  the  cycle  a 
balloon  would  follow.  (Rand,  1992,  p.  4) 

Assume  a  zero  pressure  balloon  is  floating  at  an 
equilibrium  altitude.  During  the  day,  solar  radiation  will 
heat  the  internal  gas,  resulting  either  in  increased  volume 
(if  the  balloon  is  not  fully  inflated)  or  in  decreased  density 
(if  the  balloon  is  fully  inflated,  gas  is  vented,  or  forced 
out  of  the  balloon)  .  In  either  case,  the  balloon  is  more 
buoyant  and  it  will  rise  to  a  new  density-equilibrium 
altitude.  At  night,  or  in  periods  of  cloudiness,  the  lifting 
gas  is  cooled  and  it  becomes  more  dense.  Pressure  remains 
constant,  so  the  balloon  envelope  decreases  in  volume.  The 
decrease  in  displaced  air  means  the  balloon  will  lose  altitude 
until  a  new  density-equilibrium  altitude  is  reached.  (Rand, 
1992,  p.  4) 

A  zero  pressure  balloon  may  remain  at  a  constant 
altitude  only  when  gas  is  vented  (if  the  gas  temperature  is 
rising),  or  when  ballast  is  dropped  (if  the  gas  temperature  is 
falling) .  Even  when  using  these  techniques  to  control  zero 
pressure  balloon  flight,  mission  duration  of  zero  pressure 
balloons  is  limited  to  five  to  seven  days,  except  in  rare 
environmental  conditions  such  as  at  the  poles,  because 
approximately  eight  percent  of  the  system  mass  (balloon  plus 
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fixed,  but  varies  in  altitude  depending  upon  latitude  and  time 


*t 


payload  plus  ballast)  must  be  dropped  each  night  to  maintain 

altitude.  (Rand,  1992  p.  5) 

Recent  developments  in  automatic  ballasting  have  made 

it  possible  to  extend  zero  pressure  balloon  mission  duration 

at  high  altitudes  in  unmanned  flight. 

A  group  of  researchers  has  reported  the  successful 
flight  of  a  zero  pressure  balloon  for  40  days  above 
Antarctica.  A  pressure  sensor  was  used  to 
automatically  release  two  kilograms  of  ballast 
whenever  the  balloon  descended  to  two  kilometers  below 
the  desired  float  altitude  of  30  kilometers.  (Rand, 
p.  8) 

This  example  illustrates  that  automatic  ballasting  may 

increase  mission  duration,  but  the  unique  environment  of 

Antarctica  was  a  significant  factor  in  this  case.  Thus,  it 

remains  to  be  seen  if  automatic  ballasting  would  increase 

mission  duration  so  dramatically  at  mid-latitudes. 

Zero  pressure  balloons  can  carry  several  thousand 

pounds  of  payload  to  altitudes  as  high  as  130,000  feet  because 

there  is  very  little  stress  on  the  balloon  envelope. 

In  December  1990,  a  29.47  million  cubic  foot  helium 
balloon  carried  a  3750  pound  payload  to  an  altitude  of 
130,000  over  Antarctica.  During  its  nine  day  mission, 
the  balloon  circumnavigated  the  Antarctic  continent 
(more  than  4,000  miles)  and  landed  its  payload  only 
113  miles  from  the  launch  point.  (Winzen,  1991,  p.  2) 

The  condition  of  24  hours  per  day  of  sunlight  minimized 

ballast  requirements,  gas  venting,  and  altitude  changes,  which 

maximized  the  mission  duration  of  the  balloon. 
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In  sum,  zero  pressure  balloons  can  be  used  to  reliably 
lift  large  payloads  to  significantly  high  altitudes,  but  only 
for  about  a  week  at  mid-latitudes. 

2.  Superpressure  Balloons 

Superpressure  balloons  are  similar  to  zero  pressure 
balloons,  except  that  the  envelope  is  sealed  at  the  bottom  to 
create  a  pressurized  envelope  so  no  gas  can  escape.  Even 
though  the  internal  pressure  varies  slightly  due  to  changes  in 
internal  gas  temperature,  these  balloons  remain  at  a  constant 
volume  because  of  the  strength  of  the  envelope.  This  constant 
envelope  volume  results  in  altitude  stability  at  a  constant 
density  altitude  where  system  weights  are  in  equilibrium  with 
the  surrounding  atmosphere.  The  high  internal  pressures 
require  a  film  envelope  that  is  thin  yet  strong,  lightweight, 
free  of  pinholes,  and  impermeable  to  gas  diffusion.  (Winzen, 
1991,  p.  2) 

The  first  superpressure  balloon  was  flown  shortly 
after  Montgolfier's  first  flight,  but  significant  research  of 
superpressure  balloons  was  not  pursued  until  polyethylene,  a 
synthetic  material  with  the  physical  properties  necessary  for 
superpressure  envelopes,  was  developed  during  World  War  II. 
Polyethylene  was  tested  extensively  as  the  envelope  in  zero 
pressure  balloons  during  the  1950's,  then  tried  with 
superpressure  balloons  beginning  in  1961.  From  1968  through 
1970,  numerous  superpressure  balloon  flights  were  conducted  by 
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the  National  Center  for  Atmospheric  Research  (NCAR) .  Over  200 
balloons  were  flown  at  altitudes  of  16-24  kilometers  for 
durations  of  up  to  744  days,  with  altitude  deviations  of  less 
than  100  meters.  Payloads  were  less  than  one  pound.  (Lawrence 
Livermore,  1990,  p.  52) 

In  1973,  the  National  Scientific  Balloon  Facility 

(NSBF)  conducted  Project  Boomerang: 

Two  20  meter  diameter  polyester  superpressure  spheres 
were  successfully  flown  from  Australia,  with  payloads 
of  52  kilograms  and  46  kilograms  respectively.  Both 
flights  circumnavigated  the  globe:  one  flight  lasted 
36  days  and  was  recovered  within  16  kilometers  of  the 
launch  point;  the  other  lasted  212  days,  [recovery 
status  not  stated]  (Rand,  1992,  p.  9) 

Multiple  superpressure  balloon  flights  were  conducted  during 

the  1970' s  by  various  organizations  with  mixed  results.  As 

attempts  were  made  to  scale  the  balloons  to  larger  sizes,  it 

was  found  that  the  materials  being  used  had  the  propensity  for 

catastrophic  failure  from  flaws  that  developed  either  in 

manufacture  or  in  handling.  Thus,  superpressure  ballooning 

research  was  reduced  while  zero  pressure  balloon  research  was 

emphasized.  (Winzen,  1991,  p.  1) 

Throughout  the  1970 's  and  1980' s.  Mylar  was  the  best 

material  available  for  superpressure  balloons. .. indeed  the 

small  Mylar  balloons  that  may  be  bought  at  any  number  of  gift 

or  party  stores  today  are  really  superpressure  balloons. 

Unfortunately,  Mylar  develops  pinholes  during  manufacture,  so 

it  is  a  poor  gas  barrier,  and  Mylar  superpressure  balloons  are 


limited  to  payloads  less  than  100  pounds  at  high  altitudes. 
(Winzen,  1991,  p.  1) 


Recognizing  the  potential  applications  of  stable  high 
altitude  balloons,  the  Defense  Advanced  Research  Project 
Agency  (DARPA)  sponsored  a  Small  Business  Innovative  Research 
(SBIR)  program  to  develop  a  superpressure  system  with  better 
materials.  The  best  material,  determined  from  a  1991  study  by 
Winzen  International  of  San  Antonio,  Texas,  is  a  biaxially 
oriented  nylon  film  which  has  the  strength,  weight,  and  gas 
impermeability  properties  desired.  Synthetic  films  are 
usually  stressed  in  only  one  direction  during  manufacture  but 
this  material  is  a  nylon  film  that  is  stressed  in  two 
directions.  This  process  orients  the  film's  molecules  so  that 
the  material  is  as  strong  laterally  as  it  is  longitudinally. 
Research  has  shown  that  the  altitude  of  a  balloon  made  of  this 
material  will  vary  from  day  to  night  by  only  200  meters, 
independent  of  the  size  of  the  balloon.  Furthermore, 
permeability  measurements  demonstrate  that  the  life  of  this 
balloon  should  exceed  four  years.  (Rand,  1991,  pp.  1-1  -  1-3) 
Two  test  flights  have  been  conducted  using  this 
biaxially  oriented  nylon  material.  In  August  1992,  a  nine 
meter  diameter  balloon  was  launched  in  Utah,  carrying  a  14 
kilogram  payload  to  an  altitude  of  20  kilometers.  In  October 
1992,  a  23  meter  diameter  balloon  carried  a  23  kilogram 
payload  to  an  altitude  of  33  kilometers.  Both  balloons 

maintained  their  design  altitudes  until  they  were  destroyed  by 
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command.  Test  flights  of  superpressure  balloons  with  this 
material  continue  with  the  following  goals: 


•  Within  one  year;  fly  a  23  kilogram  payload  to  36 
kilometers  for  30  days. 

•  Within  18  months;  fly  a  450  kilogram  payload  to  36 

kilometers  for  30  days.  (D.  Brown,  1992,  p.  56)  • 

3.  Sky  Anchor  Balloons 

The  sky  anchor  is  a  hybrid  system  combining  zero 
pressure  balloons  and  superpressure  balloons  in  an  attempt  to  • 

stabilize  zero  pressure  altitude  excursions  and  to  achieve 
extended  flight.  The  idea  is  to  fly  two  balloons  together  to 
gain  lift  capacity  with  the  zero  pressure  balloon  and  to  gain  • 

altitude  stability  with  the  superpressure  balloon  by  using  it 
as  air  ballast.  As  the  zero  pressure  balloon  ascends  due  to 
warming,  the  superpressure  balloon  becomes  heavier  than  the  • 

surrounding  air,  preventing  the  entire  system  from  ascending 
to  an  altitude  that  would  require  gas  to  be  vented  from  the 
zero  pressure  balloon.  Cooling  the  gas  in  the  zero  pressure  • 

balloon  returns  the  system  to  its  original  equilibrium  point. 

(Lawrence  Livermore,  1990,  p.  9) 

Sky  anchor  systems  have  been  constructed  and  flown,  • 

but  with  little  success.  The  challenges  involved  in  handling 
and  launching  two  balloons  simultaneously  are  significant,  and 
even  various  configurations  of  balloons  and  payloads  have  • 

produced  limited  results.  The  theory  of  operation  at  altitude 
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is  fine;  the  difficulty  lies  in  getting  the  system  to  an 
equilibrium  altitude.  (Winzen,  1991,  p.  3) 

The  NSBF  conducted  a  series  of  tests  with  sky  anchors 
in  the  late  1970' s.  Numerous  launch  problems  were 
experienced,  but  one  system  carrying  227  kilograms  was  able  to 
remain  at  about  36  kilometers  altitude  for  four  days. 
Unfortunately,  altitude  variations  of  up  to  six  kilometers 
were  frequent.  (Lawrence  Livermore,  1990,  p.  52) 

The  most  famous  sky  anchor  system  is  the  EarthWinds 
project,  which  in  recent  years  has  made  repeated  attempts  to 
circumnavigate  the  globe.  The  crew's  attempts  have  not  been 
successful,  and  they  have  encountered  significant  skepticism 
from  the  scientific  community  and  the  media  throughout. 
Overall,  the  sky  anchor's  poor  record  make  it  an  unlikely 
candidate  for  military  use.  (S.  Brown,  1992,  pp.  80-126) 

C.  SUMMARY  OF  BALLOON  USE  AMD  HISTORY 

A  variety  of  balloons  have  been  used  effectively  for  a 
number  of  purposes  in  military  conflicts  since  their  invention 
by  Frenchmen  in  1782.  Two  types  of  balloons,  zero  pressure 
and  superpressure,  seem  to  offer  good  potential  for  future 
military  use.  Currently,  zero  pressure  balloons  are  capable 
of  carrying  large  payloads  for  up  to  a  week  while 
superpressure  balloons  are  capable  of  smaller  payloads  for 
many  weeks  or  even  years. 


•  • 
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IV.  ATMOSPHERIC  CIRCOLATIOM  AMD  MODELING 

Central  to  any  discussion  of  the  use  of  free-floating 
balloons  is  an  understanding  of  atmospheric  dynamics.  This 
section  describes  the  differences  in  the  earth's  atmospheric 
circulation  based  upon  altitude,  latitude,  and  time  of  year. 


A.  ATMOSPHERE  STROCTORE 

Meteorologists  conventionally  divide  the  atmosphere  into 
four  layers  based  on  the  vertical  gradient  of  temperature. 
These  layers,  as  shown  in  Figure  2,  are  the  troposphere, 
stratosphere,  mesosphere,  and  the  thermosphere.  Temperature 


Figure  2.  U.S.  standard 
Atmosphere  Temperature 
Profile,  1976 


generally  decreases  with  height  in  the  troposphere  and  in  the 
mesosphere,  and  it  generally  increases  with  height  in  the 
stratosphere  and  in  the  thermosphere.  The  troposphere 
accounts  for  almost  85  percent  of  the  total  mass  of  the 
atmosphere,  and  it  accounts  for  virtually  all  atmospheric 
water  vapor.  The  tropopause,  which  separates  the  troposphere 
and  the  stratosphere,  is  a  level  of  temperature  minimum  which 
varies  in  height  from  about  15  kilometers  at  the  equator  to 
nine  kilometers  at  the  poles.  The  stratopause  is  a  level  of 
temperature  maximum  near  50  kilometers  which  separates  the 
stratosphere  and  the  mesosphere.  The  mesosphere  is  bounded 
above  by  the  mesopause,  a  level  of  temperature  minimum  similar 
to  the  tropopause  at  about  80  kilometers.  (Andrews,  1987,  p.3) 
In  addition  to  temperature,  three  other  physical 
properties  characterize  the  atmosphere:  pressure,  density  and 
velocity.  Pressure  is  defined  as  the  amount  of  force  applied 
over  a  surface  or  force  per  unit  area,  which  can  be  either  the 
earth's  surface  or  an  air  parcel.  Pressure  is  measured  in 
either  millibars  (mb)  or  in  kilopascals  (kPA) ,  with  one  bar 
equal  to  14.5  pounds  per  square  inch.  Standard  sea  level 
pressure  (STP)  equals  1013.25  mb  or  101.325  kPa.  Air  pressure 
decreases  with  altitude;  several  altitude-pressure  reference 
points  are  provided  in  Table  III.  Density  is  the  amount  of  a 
substance  per  unit  measure  or  mass  at  standard  pressure  and 
temperature.  The  air  density  also  decreases  with  altitude  as 
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TABLE  III.  ALTITUDE-PRESSURE 
REFERENCE  POINTS 


both  pressure  and  temperature  decrease.  Velocity  is  a  measure 
of  wind  speed  relative  to  the  ground,  measured  in  meters  per 
second.  As  noted  in  the  discussion  on  balloon  buoyancy,  all 
three  of  these  physical  quantities  have  an  impact  upon  balloon 
lift  and  drift.  (Holton,  1979,  p.l) 


B.  ATMOSPHERIC  CIRCULATION 


1.  Fundaaental  Motion  Forces 


The  motions  of  the  atmosphere  are  governed  by  the 
fundamental  laws  of  fluid  mechanics  and  thermodynamics:  the 
laws  of  conservation  of  mass,  momentum,  and  energy.  The 
primary  forces  which  cause  atmospheric  motion  are  the  pressure 
gradient  force,  the  gravitational  force,  and  friction. 
Additionally,  one  "apparent"  force  also  acts  upon  the 
atmosphere  to  cause  motion:  the  Coriolis  force.  (Holton, 
1979,  p.5-17) 

An  object  or  air  parcel  accelerates  in  the  direction 
of  the  applied  pressure  gradient  force,  as  shown  in  Figure  3. 


The  pressure  gradient  force  is  proportional  to  the  gradient  of 
the  pressure  field,  not  to  the  pressure  itself.  This  force  is 
active  through  the  entire  atmosphere.  (Holton,  1979,  p.5-7) 
Tropospheric  and  stratospheric  circulation  are  both 
strongly  influenced  by  the  Coriolis  force,  which  is  an 
apparent  force  or  effect  which  accounts  for  the  rotation  of 
the  earth.  The  rotation  of  the  earth  imparts  spin  to  an  air 
particle  in  the  atmosphere,  which  causes  the  air  particle  to 
have  an  angular  momentum,  or  coriolis  acceleration,  with 
respect  to  the  earth.  As  shown  in  Figure  4  (after  Lawrence 


Figure  4.  Coriolis 
Acceleration 
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Livermore) ,  this  angular  momentum  varies  with  latitude,  with 
no  Coriolis  acceleration  for  horizontal  motion  at  the  equator. 
As  the  air  particle  is  accelerated  by  the  pressure  gradient 
force,  the  Coriolis  acceleration  will  increase  until  it 
reaches  a  balance  as  shown  in  Figure  5  (after  Lawrence 


Figure  5.  Geostrophic  Winds 


Livermore) .  The  balance  causes  "geostrophic"  winds  that  flow 
parallel  to  the  isobars  at  a  speed  proportional  to  the 
pressure  gradient  and  inversely  proportional  to  the  sine  of 
the  latitude.  (Lawrence  Livermore,  1990,  p.  25) 

Friction  is  important  near  the  surface  of  the  earth, 
where  atmospheric  motion  is  retarded  by  contact  with  the 
earth's  surface.  Friction  thus  influences  circulation  in  the 
lower  troposphere  but  its  only  effect  at  higher  altitudes  is 
indirect  through  the  interaction  of  tropospheric  eddy  motions 
with  the  lower  stratosphere.  (Holton,  1979,  pp.  296-298) 


4r, 


■ 


26 


•  •  • 


2.  Oba«rv«d  Circulation 

a.  Zonal ly  Averaged  Circulation 

All  references  in  atmospheric  circulation  studies 
create  a  distinction  between  the  longitudinally  averaged  flow, 
which  is  either  zonal  mean  flow  or  meridional  mean  flow,  and 
the  deviations  from  these  means,  or  eddies.  Zonal  mean  winds 
are  parallel  to  the  equator,  while  meridional  mean  winds  are 
perpendicular  to  the  equator.  This  section  describes 
stratospheric  circulation  in  terms  of  zonal ly  averaged 
circulation. 

b.  Observed  Circulation  Patterns 

A  combination  of  radiosonde  data,  rocketsonde  data, 

and,  more  recently,  remote  temperature  soundings  from 

satellites,  have  provided  meteorologists  with  a  much  clearer 

picture  oi  stratospheric  dynamics  than  ever  before.  At  this 

time,  the  study  of  upper  atmosphere  circulation  patterns  has 

revealed  that  there  are  only  a  handful  of  general  cases  of 

circulation,  and  that  upper  atmosphere  circulation  is 

influenced  by  the  annual  solar  cycle. 

The  net  radiative  heating  distribution  has  a  strong 
seasonal  dependence  with  maximum  heating  at  the  summer 
pole  and  maximum  cooling  at  the  winter  pole.  The  CSoriolis 
torque  exerted  by  this  meridional  flow  generates  mean 
zonal  easterlies  (from  the  east)  in  the  summer  hemisphere 
and  westerlies  in  the  winter  hemisphere.  (Andrews,  1987, 

p.  6) 

The  circulation  patterns  vary  gradually  from  month  to 
month  during  the  annual  cycle  and  they  recur  regularly.  The 
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general  categories  of  stratospheric  circulation  are  listed 
below: 

•  the  extratropical,  or  non-equatorial,  pattern  is  westerly 
zonal-mean  winds  in  the  winter  hemisphere,  and  easterly 
zonal-mean  winds  in  the  summer  hemisphere.  Figure  6  shows 
zonal  mean  winds  in  meters  per  second  for  solstice 
conditions  with  W  and  E  designating  centers  of  westerly 
(positive,  from  the  west)  and  easterly  (negative,  from  the 
east)  winds,  respectively.  (Andrews,  1987,  p.  8) 


* 


♦ 


•  stratospheric  sudden  winter  warmings  in  the  Northern 
Hemisphere  lead  to  mean-flow  deceleration  and  even  to  a 
reversal  of  the  winds  to  easterly.  (Andrews,  1987,  p. 
259) 

•  the  equatorial  pattern  is  an  alternating  pattern  of 
eastward  and  westward  winds  that  repeat  at  intervals 
varying  from  about  22  to  34  months,  with  an  average  period 
of  about  27  months  (though  there  is  a  six  month  cycle  at 
higher  levels).  (Andrews,  1987,  p.  313) 

•  the  transition  between  the  Northern  and  the  Southern 
Hemispheres  at  the  equinoxes  results  in  weak  mean  zonal 
westerlies  in  both  hemispheres.  (Andrews,  1987,  p.  6) 

An  important  point  to  observe  from  Figure  6  is  that 
wind  speed  seems  to  be  at  a  minimum  at  approximately  2  0 
kilometers  altitude.  This  region  of  minimum  winds  is  not 
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fixed,  but  varies  in  altitude  depending  upon  latitude  and  time 
of  year.  Figures  7-10  show  that  altitudes  of  30-40  kilometers 
appear  to  be  the  most  promising  for  balloon  surveillance 
operations. 


Figure  10.  October  Zonal  Mean  winds 


5.  Real  World  Meaning  of  Results 


C.  ATMOSPHERIC  CIRCOLATIOM  MODELING 


The  majority  of  atmospheric  circulation  models  that  have 
been  developed  are  tropospheric  models,  primarily  because  of 
the  need  to  understand  the  earth's  climate  and  weather. 
However,  much  more  work  dealing  with  the  stratosphere  has  been 
done  in  the  last  20-30  years  because  of  concerns  about  the 
earth's  ozone  layer.  Circulation  models  fall  into  two 
categories:  those  based  on  simulated  physical  processes  and 
those  based  on  historical  data. 

One  general  circulation  model  based  on  a  simulation  of  the 
physical  processes  is  the  SKYHI  model,  developed  by  the 
Geophysical  Fluid  Dynamics  Laboratory  at  Princeton.  SKYHI  is 
considered  to  provide  the  most  complete  representation  of  the 
middle  atmosphere:  it  has  40  prediction  levels  extending  from 
the  surface  of  the  earth  to  about  80  kilometers  and  the 
predicted  fields  include  wind,  temperature,  water  vapor  and 
surface  pressure.  (Andrews,  1987,  p.  422) 

A  model  based  on  historical  atmospheric  data  is  the  Global 
Reference  Atmosphere  Model,  or  GRAM,  an  empirical  FORTRAN 
computer  simulation  of  the  earth's  atmosphere  developed  at 
Georgia  Tech.  GRAM,  which  is  currently  the  responsibility  of 
the  National  Aeronautical  and  Space  Administration  (NASA) , 
uses  data  from  several  sources  and  over  several  years  to 
compute  monthly  averages  of  atmospheric  properties,  including 
wind  speed  and  wind  direction.  GRAM  provides  a  worldwide,  12- 


month  database  of  properties  in  four  dimensions:  latitude, 
longitude,  altitude,  and  time  of  year.  There  are  two  key 
points  to  keep  in  mind  when  using  the  GRAM.  The  first  is  that 
the  GRAM  contains  statistically  averaged  data:  it  provides 
mean  wind  speed  and  direction,  with  standard  deviations,  for 
an  entire  month  for  a  given  combination  of  latitude, 
longitude,  and  altitude.  Secondly,  insufficient  data  for  the 
Southern  Hemisphere  caused  the  designers  of  the  GRAM  to 
represent  it  with  a  six  month  displacement  of  Northern 
Hemisphere  data.  (Hawkins,  1991,  pp.  1-18) 

This  second  shortcoming  of  GRAM  was  alleviated  in  its 
latest  version,  produced  in  1990  and  appropriately  called 
GRAM-90.  GRAM-90  incorporates  extensive  new  data,  mostly 
gathered  by  satellites,  to  utilize  actual  data  from  the 
Southern  Hemisphere  for  each  month.  Like  GRAM,  GRAM-90 
provides  mean  values,  but  it  also  produces  "perturbation" 
values  which  are  slightly  different  from  the  mean  values. 
These  perturbation  values  are  drawn  from  the  historical  record 
of  winds  to  provide  a  broad  perspective  of  atmospheric 
conditions  at  a  given  point.  (Jeffries,  1993) 

Other  models  are  available,  but  the  models  just 
described  demonstrate  the  basic  types  and  capabilities  of 
circulation  models.  Of  these,  SKYHI  is  more  of  a  "prediction" 
model,  while  GRAM,  and  the  GRAM-90  in  particular,  provide 
useful  databases  for  statistical  analysis. 
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narrow  over  short  periods  of  time,  a  distribution  analysis 


D.  80NNARY  OF  ATMOSPHERIC  CXRCOLATXOM  AMD  MODELING 


To  sum  up,  it  should  be  noted  that  stratospheric  • 

circulation  is  different  than  that  of  the  troposphere  and  that 
models  which  describe  atmospheric  circulation  are  available. 
Stratospheric  winds  are  fairly  stable,  relatively  light  and  • 

their  circulation  patterns  vary  during  the  year  as  a  result  of 
the  annual  solar  cycle.  However,  even  though  wind  circulation 
seems  favorable  for  military  balloon  operations,  it  remains  to  • 

be  seen  if  the  circulation  patterns  will  cause  a  free  floating 
balloon  to  "orbit”  the  earth  or  if  they  will  cause  them  to 
congregate  in  localized  patterns.  This  is  the  subject  of  the  • 

next  chapter. 
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V.  TRAJECTORY  PREDICTION 


Generally  speaking,  people  reject  out  of  hand  the  idea  of 
using  free-floating  balloons  for  military  surveillance  or 
communications  missions  because  of  the  uncertainty  involved. 
Balloon  trajectories  cannot  be  predicted  accurately  like 
satellite  orbits,  nor  are  balloons  tightly  controlled  like 
aircraft.  However,  the  variability  of  free-floating  high 
altitude  balloon  trajectories  is  not  extreme:  as  described  in 
Chapter  IV,  winds  in  the  upper  stratosphere  may  be  stable 
enough  to  provide  fairly  constant  trajectories.  Balloons 
travelling  in  these  winds  may  even  be  able  to  provide  better 
coverage  of  an  area  than  low-earth  orbiting  satellites  because 
of  longer  time  over  target.  Thus,  a  statistical  analysis  of 
the  variability  of  balloon  trajectories  in  the  upper 
stratosphere  may  be  useful  either  to  confirm  or  to  refute 
conventional  wisdom  about  the  potential  of  free-floating 
balloons. 

A.  BALLOON  TRAJECTORY  PREDICTION  PROGRAMS 

A  study  of  relevant  literature  reveals  that  very  little 
work  has  been  done  in  the  area  of  balloon  trajectory 
prediction,  with  the  exception  of  predicting  vertical 
trajectories  of  weather  balloons.  The  paucity  of  study  led 
DARPA  to  sponsor  Coleman  Research  Corporation  (CRC)  to  develop 
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a  Balloon  Drift  Pattern  Simulation  (BDPS) .  DARPA  personnel 
were  interested  in  exploiting  high-altitude,  expendable 
balloons  for  a  variety  of  military  missions,  and  they  needed 
a  model  to  predict  how  balloons  would  drift  in  order  to  assess 
the  feasibility  of  several  concepts. 

In  Phase  I  of  an  SBIR  project,  CRC  wrote  an  upper- 
atmosphere  drift  pattern  simulation  for  execution  and  display 
output  on  VAX  computers,  demonstrating  the  technical 
feasibility  of  predicting  drift  patterns  using  a  digital 
computer  simulation.  In  Phase  II  of  the  SBIR,  DARPA  required 
CRC  to  develop  a  Macintosh-based  version  of  the  BDPS.  The 
intent  was  to  build  a  desk-top,  deployable  aid  that  would 
enable  a  theater-level  commander  to  assess  how  to  employ 
balloon  systems  for  either  communications  or  surveillance 
missions.  (Hawkins,  1991,  p.  1-18) 

The  BDPS  is  a  time-step  simulation  which  draws  upon  a 
database  of  wind  tables  to  compute  a  balloon  drift  pattern. 
This  database  can  be  actual  wind  tables  (either  archived  wind 
data  or  forecast  wind  data)  or  a  climate  model.  The  user 
chooses  the  database  based  on  his  needs:  archived  and  forecast 
data  (such  as  the  Navy's  NOGAPS  data  from  the  Fleet  Numerical 
Oceanography  Center)  is  only  available  to  an  altitude  of  about 
30  kilometers  (10  millibars).  If  the  simulation  calls  for 
drift  patterns  above  this  altitude,  the  forecast  data  must  be 
extrapolated  or  a  climate  model  must  be  used.  As  described  in 
Chapter  IV,  several  climate  models  are  available  for  use. 


•  • 
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since  BDPS  is  used  on  a  Macintosh,  it  is  mostly  user-friendly, 
but  it  is  quite  large  (39  Megabytes)  and  it  is  also  very  slow: 
it  takes  about  10  minutes  to  calculate  and  to  display  just  24 
hours  of  a  trajectory.  (Hawkins,  1991,  p.  1-18) 

The  combination  of  the  BDPS  program  and  the  DARPA  study  on 
superpressure  balloons  demonstrated  the  potential  of  balloon 
systems.  Unfortunately,  the  personnel  at  DARPA  who  were 
interested  in  such  systems  rotated  to  other  duty  stations  and 
the  balloon  system  ideas  have  not  been  pursued. 

The  GRAM-90  is  not  a  drift  prediction  model  in  that  it 
provides  atmospheric  data  for  a  specific  point  of  altitude, 
latitude,  longitude,  and  date.  However,  modifications  were 
made  to  the  program  so  that  it  could  be  used  as  to  describe 
balloon  drift.  Specifically,  source  code  was  added  to  convert 
the  wind  direction  and  wind  speed  for  a  designated  time  step 
into  a  drift  vector  to  apply  to  the  initial  balloon  start 
point  (altitude,  latitude,  longitude,  and  date) .  This  vector 
was  used  to  compute  a  new  balloon  location.  This  new  location 
was  then  treated  as  another  start  point  input  for  the  model, 
which  produced  the  atmospheric  conditions  at  the  new  location. 
This  wind  vector  and  time  step  process  was  repeated  for  a  year 
and  the  balloon's  location  was  recorded  intermittently  for 
future  analysis.  The  design,  use  and  analysis  of  this 
simulation  are  described  in  more  detail  in  the  next  section. 
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B.  PREDICTION  SIMULATION  DESIGN 

The  basic  design  of  the  simulation  was  to  use  the  modified 
GRAM-90  as  a  time-step  simulation  of  balloon  movement 
resulting  from  winds  at  altitude.  By  changing  the  input  seed 
to  a  random  number  generator  used  in  the  unmodified  program, 
varying  trajectories  were  produced.  The  seed  was  randomly 
changed  for  each  replication  of  a  given  set  of  initial 
conditions  (latitude,  longitude,  altitude  and  month) .  Balloon 
locations  were  recorded  as  desired  to  provide  data  for  a 
determination  of  the  variability  of  the  trajectories. 

1.  Physical  Equipment 

The  modified  GRAM-90  was  used  on  the  Naval 
Postgraduate  School's  Amdahl  mainframe  computer  emulating  an 
IBM  3270.  Data  analysis  was  completed  using  the  Minitab 
statistical  package,  also  on  the  mainframe. 

2 .  Bypotheses 

The  null  hypothesis  was  that  the  trajectories  of  free- 
floating  balloons  drifting  at  36  kilometers  altitude  would  not 
follow  regular  or  repeating  patterns  of  drift.  The 
alternative  hypothesis  was  that  these  balloons  would  indeed 
follow  regular  or  repeating  patterns  of  drift. 

3.  Assumptions 

•  GRAM-90  as  modified  is  a  valid  balloon  drift  prediction 
simulation.  Since  a  balloon  will  achieve  equilibrium  with 
the  winds  surrounding  it,  the  assumption  is  that  a  one 
hour  time  step  is  sufficiently  small  to  accurately 
describe  the  balloon's  trajectory  using  only  the  wind 
speed  and  wind  direction  outputs  from  the  model. 
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•  GRAM-90  is  a  valid  atmospheric  model  of  winds  at  altitude. 
Although  the  GRAM-90  has  averaged  data,  the  assumption  is 
that  drawing  the  standard  deviations,  or  perturbation 
values,  from  many  years  of  data  ensured  validity.  In 
fact,  the  wind  speeds  varied  considerably  more  than 
expected,  although  they  tended  to  average  out  over  time. 

•  Balloon  altitude  is  constant.  The  vertical  position  of 
the  balloon  was  assumed  to  be  unchanging  in  an  attempt  to 
study  the  horizontal  trajectory  of  a  drifting  balloon. 

•  Balloon  ascent  is  assumed  to  be  instantaneous  to  a  point 
directly  overhead  the  starting  latitude  and  starting 
longitude.  GRAM-90  describes  atmospheric  conditions  at  a 
specific  point;  the  available  portion  was  conditions 
above  30  kilometers.  Tropospheric  winds  which  would  have 
an  effect  on  the  horizontal  position  of  a  balloon  during 
ascent  were  not  be  modeled  for  this  simulation. 

4.  Measure  of  Balloon  Drift 

Balloon  location  was  measured  in  latitude  and 
longitude  after  ten  days,  30  days,  and  360  days.  Latitude  was 
in  a  form  suitable  for  analysis  in  its  range  from  -90°  to  +90°. 
Longitude  was  not  in  a  suitable  form:  in  normal  mathematical 
operations  -179  is  178  units  away  from  +179,  yet  there  is  a 
difference  of  only  two  degrees  between  -179°  and  +179°.  Thus, 
it  was  necessary  to  convert  longitude  from  absolute  degrees  to 
a  degree  difference  from  the  starting  longitude. 

5.  Statistical  Design  of  Simulation 

Table  IV  lists  the  simulation's  variables  and  values, 
which  were  chosen  since  atmospheric  conditions  vary  with  each. 


TABLE  IV.  SIMULATION  VARIABLES 
AND  VALUES 


Month 

Latitude 

Lonoitude 

March 

45°  S 

90°  E 

June 

0° 

90°  W 

September 

45°  N 

December 

•  •  • 
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The  values  of  each  variable  are  not  all-inclusive  but  were 
chosen  to  provide  a  cross-section  for  each  factor.  Each  of 
the  24  possible  combinations  of  these  variables  was  used  to 
specify  the  initial  conditions  for  multiple  replications  of 
the  simulation. 

C.  PREDICTION  SIMULATION  DATA  ANALYSIS 

Data  was  separated  into  three  categories  for  analysis:  ten 
days;  30  days;  and  360  days.  Data  in  each  category  was 
checked  for  normality  and  uniformity,  and  then  the  variance  of 
the  means  and  medians  was  compared  through  a  variety  of 
statistical  tests.  Lastly,  the  three  categories  were  compared 
to  each  other  for  a  broader  perspective. 

1.  Balloon  Location  After  Ten  Days 

Data  from  all  24  sets  of  initial  conditions  in  this 
category  was  found  to  be  neither  normal  nor  uniform.  Latitude 
of  a  balloon  after  drifting  ten  days  does  follow  a  pattern  but 
longitude  does  not.  Despite  a  lack  of  normality,  an  Analysis 
of  Variance  (ANOVA)  for  latitude  shown  in  Table  V  was  used  to 
indicate  which  factors  were  significant.  Since  the  data  was 
not  normal,  the  numerical  values  derived  from  the  ANOVA  were 
not  valid.  The  ANOVA  shows  that  ending  latitude  was  dependent 
upon  starting  latitude  and  month,  but  not  on  starting 
longitude  (P  =  0.547).  Ending  longitude  was  also  dependent 
upon  starting  latitude  and  month,  and  it  was  dependent  upon 
starting  longitude  only  for  Equator  launches. 
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LATITUDE  ANALYSIS  OF  VARIANCE  AT  TEN  DAYS 


MTB  >  anova  ' lat 

'  ”  'month' 1 'fllat*  J 

'slong' 

Factor 

Type 

Levels  Values 

MONTH 

fixed 

4 

3  6 

9 

12 

SLAT 

fixed 

3 

-45  0 

45 

SLONG 

fixed 

2 

-90  90 

Analysifl 

of  Variance  for 

LAT 

Source 

DF 

SS 

MS 

F 

P 

MONTH 

3 

50405 

16802 

521.40 

0.000 

SLAT 

2 

1158260 

579130 

1.8E*04 

0.000 

SLONG 

1 

12 

12 

0.36 

0.547 

MONTH* SLAT 

6 

15835 

2639 

81.90 

0.000 

MONTH*SLONG 

3 

5394 

1798 

55.79 

0.000 

SLAT •SLONG 

2 

4755 

2377 

73.77 

0.000 

MONTH  *  SLAT  *  SLONG 

6 

5796 

966 

29.98 

0.000 

Error 

720 

23202 

32 

Total 

743 

1263658 

LAT 

ENDING  LATITUDE 

DF  «  Degrees  of 

Freedom 

MONTH  » 

START  MONTH 

SS  >  Sum  of  Squares 

SLAT  - 

START  LATITUDE 

MS  >  Mean  Square 

SLONG  « 

START  LONGITUDE 

F  =  F 

Test  Statistic 

P  =  P- 

Value 

Since  longitude  was  found  not  to  be  a  significant  factor,  data 
was  analyzed  further  in  the  configuration  shown  in  Table  VI. 


TABLE  VI . 


INITIAL  CONDITION  SETS 


SET  NUMBER  MONTH  LATITUDE  LONGITUDE 

1 

MAR 

45°  S 

90°  E/90°  W 

2 

MAR 

0° 

90°  E/90°  W 

3 

MAR 

45°  N 

90°  E/90°  W 

4 

JUN 

4  5°  S 

90°  E/90°  W 

5 

JUN 

0° 

90°  E/90°  W 

6 

JUN 

45°  N 

90°  E/90°  W 

7 

SEP 

45°  S 

90°  E/90°  W 

8 

SEP 

0° 

90°  E/90°  W 

9 

SEP 

45°  N 

90°  E/90°  W 

10 

DEC 

45°  S 

90°  E/90°  W 

11 

DEC 

0° 

90°  E/90°  W 

12 

DEC 

45°  N 

90°  E/90°  W 
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Figure  11  shows  mean  latitudinal  displacement  distances  for 
this  category  where  balloons  remained  within  about  eight 
degrees  of  their  starting  latitude.  This  deviation  varies  by 


Figure  11.  10  Day  Latitude 

Displacement 


month,  with  June  (initial  condition  sets  4-6)  being  the  month 
when  the  balloons  remained  closest  to  the  starting  latitude 
and  March  (sets  1-3)  having  the  largest  deviation. 

2.  Balloon  Location  After  30  Days 

As  in  the  ten  day  category,  data  was  found  to  be 
neither  normal  nor  uniform.  Figure  12  shows  mean  latitudinal 


Figure  12.  30  Day  Latitude 

Displacement 


displacement  distances  for  this  category,  where  balloons 
drifted  within  about  nineteen  degrees  of  their  starting 
latitude.  The  deviations  were  smallest  again  in  June.  ANOVA 
indicated  that  balloon  location  after  30  days  is  dependent 
upon  all  factors:  month,  starting  latitude,  and  starting 
longitude. 

3.  Balloon  Location  After  360  Days 

Figure  13  shows  mean  latitudinal  displacement  for 
balloons  after  360  days  of  drift:  the  balloons  usually  drifted 


Figure  13.  360  Day  Latitude 
Displacement 

over  sixty  degrees  from  their  starting  latitude.  In  fact. 
Table  VII  shows  that  balloons  tend  to  drift  toward  the  poles 
after  almost  a  year.  Once  again,  the  data  is  neither  normal 
nor  uniform.  The  ANOVA  for  latitude  again  indicates  that 
ending  latitude  is  not  related  to  starting  longitude,  but  it 
is  dependent  upon  both  starting  latitude  and  the  starting 
month.  Most  often,  balloons  launched  from  North  of  the  Equator 
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TABLE  VIZ.  360  DAY  LATITUDE 
DISTRIBUTION 


Histogram  of 

LAT 

N  •  744 

Each  *  rtpresants 

10  obs. 

Midpoint  Count 

-80 

350 

-60 

9 

* 

-60 

8 

* 

-20 

36 

m*m* 

0 

26 

20 

38 

40 

19 

•« 

60 

88 

80 

170 

tended  to  drift  toward  the  North  Pole,  and  balloons  launched 
from  South  of  the  Equator  tended  to  drift  toward  the  South 
Pole,  although  there  were  anomalies  where  balloons  drifted 
toward  the  pole  opposite  their  launch  hemisphere. 

4.  Comparison  of  Categories 

It  can  be  seen  from  all  three  categories  that  balloons 
follow  regular  or  repeating  patterns  in  their  latitudinal 
displacement.  Further,  as  might  be  expected,  the  displacement 
from  start  latitude  increased  over  time.  In  all  three 
categories,  ending  latitude  was  found  to  be  dependent  upon 
starting  latitude  and  month  but  not  on  starting  longitude. 
Starting  longitude  was  found  not  to  be  a  significant  factor, 
except  for  an  Equator  launch  in  the  30  day  category. 


5.  Real  World  Meaning  of  Results 

The  most  important  finding  of  the  analysis  of  the  data 
is  that  balloon  drift  patterns  over  ten  days  follow  fairly 
narrow  patterns.  An  examination  of  the  distances  involved 
illustrates  this  importance.  After  ten  days,  most  balloons 
completed  between  one-half  global  circumnavigation  and  one 
complete  global  circumnavigation.  The  mean  difference  in 
latitude  at  that  time  was  8.5  degrees,  or  a  ground  distance  of 
approximately  950  kilometers.  At  an  altitude  of  36 
kilometers,  the  balloon  has  a  line-of-sight  footprint  on  the 
ground  of  730  kilometers  radius  or  1460  kilometers  diameter. 
Thus,  even  with  a  drift  distance  of  950  kilometers,  the 
balloon  retains  visibility  of  a  sizable  portion  of  the 
original  footprint,  although  that  footprint  obviously  would 
have  shifted  dramatically  longitudinally.  A  "theater" 
surveillance  balloon  which  would  be  employed  for  less  than  ten 
days,  perhaps  even  for  as  little  as  one  or  two  days,  would 
have  even  less  drift. 


VI .  CONCLOSIOMS 


A.  BALLOON  SURVEILLANCE  SYSTEM  FEASIBILITY 

The  narrow  drift  patterns  of  balloon  flights  of  ten  days 
support  the  idea  of  "theater"  or  "tactical"  high  altitude 
surveillance  balloons,  probably  with  flights  of  three  days  or 
less.  Such  balloon  systems  could  be  recovered  by  steerable 
parachute  or  mid-air  "snatch"  after  theater  transit.  With  a 
ground  speed  much  less  than  that  of  satellites  and  if 
outfitted  with  appropriate  sensors,  they  would  provide  a 
surveillance  system  that  could  overfly  enemy  areas  to  identify 
and  to  locate  mobile  ballistic  missile  systems  to  be  attacked. 
Such  a  system  would  provide  a  significant  improvement  in  the 
United  States'  missile  defense  posture. 


B.  AREAS  FOR  FUTURE  STUDY 

There  are  many  potential  areas  for  further  study  with  this 
subject.  One  suggestion  would  be  to  actually  launch  balloons 
to  check  the  validity  of  the  simulation  model.  Also,  other 
simulations  would  be  valuable.  It  is  recommended  that  the 
simulation  be  replicated  at  different  altitudes.  This  study 
looked  only  at  36  kilometers;  it  would  be  helpful  to  examine 
balloon  drift  at  a  variety  of  altitudes  to  examine  the 
variance  based  on  altitude.  Additionally,  and  perhaps  more 
importantly,  since  it  appears  that  drift  patterns  are  fairly 
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narrow  over  short  periods  of  time,  a  distribution  analysis 
similar  to  this  study  should  be  completed  with  balloon 
locations  recorded  at  intervals  less  than  ten  days.  Another 
area  to  be  looked  at  would  be  a  time  series  analysis  of 
balloon  location  to  try  to  get  more  of  a  "continuous" 
perspective  rather  than  the  discrete  approach  used  here.  Many 
other  subjects  for  study  fall  under  the  heading  of  operational 
requirements  of  a  high  altitude  balloon  system:  sensor 
selection,  payload  configuration,  power  requirements  and 
sources,  C3I  architecture,  and  concept  of  operations. 


APPENDIX  (MODIFIED  6RAM-90  800RCE  CODE) 


The  Scientific  Model  (SCIMOD  or  SCIM)  is  only  a  fraction 
of  GRAM-90,  but  since  it  is  the  only  section  in  which 
modifications  were  made  for  the  balloon  trajectory,  it  is  the 
only  portion  enclosed  for  future  reference.  The  entire  GRAM- 
90  may  be  obtained  from  NASA  Marshall  Space  Flight  Center. 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


SUBROUTINE  SCIHOO(NPOP) 

..COMPUTES  VALUES  P,D,T,U,V  AND  SHEAR  DUH.OVH  FROM  INPUT  AND 
ARRAYS  IN  COMMON  POTCOM.  INPUT  TO  SCIMOO  IS 
G  >  GRAVITY  AT  POSITION  RI  >  RADIUS  AT  HEIGHT  H 

PHIR  >  LATITUDE  (RADIANS)  THETR  •  LONGITUDE  (RADIANS) 
F10  >  F10.7  SOLAR  FLUX  F10B  «  MEAN  F10.7  FLUX 

AP  •  SOLAR-GEOMAGNETIC  A  SUB  P  INDEX 
MN/IOA/IYR  >  DATA  (lYR  <  FULL  YEAR-1900) 

IHR  MIN  <  TIME  HI  >  PREVIOUS  HEIGHT 

PHI1R  <  PREVIOUS  LATITUDE  THET1R  >  PREVIOUS  LONGITUDE 
RP1,R01,RT1  s  PREVIOUS  RANDOM  PERTURBATIONS 
SP1,S0I,ST1  *  PREVIOUS  RANOOM  STANOARD  DEVIATIONS  (SIGMAS) 
RUI.RVI  •  PREVIOUS  RANOOM  WINDS 
SUl.SVI  <  PREVIOUS  RANOOM  UINO  SIGMAS 
COMMON/IPRTP/  IPRT.NLIMIT 

C0MMON/IOTEHP/lOTEH1,I0TEM2,IUS,00,XMJO,PHI1,PHI, 

-NSAME.RPU, ROIL, RTU.SPU, SOIL, ST1L,RU1L,RV1L,SU1L,SV1L, 

S  MN,  IDA,  lYR,  HI,  PHI1R,THET1R,G,RI,H, PHIR, THETR, Ft0,F10B,AP, 
IHR,MIN,NMORE,DX,HL,VL,OZ,B, EPS, lOPP, LOOK, lET, FLAT, 
1RP1S,R01S,RT1S,RU1S,RV1S,SP1S,S01S,ST1S,SU1S,SV1S, 

2UDS 1 , V0S1 , UDL 1 , VDL 1 , UDS2 , VDS2 , U0L2 , V0L2 
COMMON  /POTCOM/IUG, MONTH, lOPR, 

.  PSP(15,19,18),DSP(15,19,18),ISP(15,19,18),USP(15,19,18>, 

.  VSP(15,19,18), 

.  PG(21,19),0G(21,19),TC(21,19),UG(21,19), 

.  PA0(17,5),0A0(17,5),TAQ(17,5),UA0(17,5),VAQ(17,5), 

.  P00(17,5),000(17,5),TDQ(17,5),U00(17,5),VDO(17,5), 

.  PH(29,19),0R(29,19),TR(29,19),UR(29,19),VR(29,19) 

.  ,Pa,OQ,TQ,UO,VQ,PaA,DOA,TaA,UA,VA,IOPa, 
1PLP(25,10),OLP(25,10),TLP(25,10), 
2ULP(25,10),VLP(25,10),U0L(25,10). 
3VDL(2S,10),U0S(2S,10),VDS(2S,10) 

COMMON  /C4/  GLAT(16),GL0N(16),NG,PGD(16, 26), 040(16,26), T4D(16, 
.  SP4(16,26)  S04(16,26),ST4(16,26),TNET1,THET,DUHMY 
COMMON/COMPEN/SPH,SOH,STH,PRH,DRH,TRH,URH,VRH,SUH,SVH,CP, 
1PRHS,DRHS,TRHS,URHS,VRHS,PRHL,DRHL,TRHL,URHL,VRHL, 
2SPHS,S0HS,STHS,SUHS,SVNS,SPHL,S0HL,STHL.SUHL,SVHL 
COMMON/UI NCOM/DH ,  FCORY  ,0X5  ,DYS  ,DPX  ,0PY  ,UGH,  VGH, 

S  TH,OTX,OTY,DUK,DVH,PH,UPRE,VPRE,OUPRE,OVPRE 
COMMON/CHK/PCK(4,4,3),OCK(4,4,3),NO(2) 

COMMON  /CHIC/LA(4,4),NB(2),IUSYM,USH,VSH,DUSH,DVSH 
COMMON  /VERT/RU1,SW1,URH,SUH,UR(29) 


SCIM 

1 

SCIM 

2 

SCIM 

3 

SCIM 

4 

SCIM 

5 

SCIM 

6 

SCIM 

7 

SCIM 

8 

SCIM 

9 

SCIM 

10 

SCIM 

11 

SCIM 

12 

SCIM 

13 

SCIM 

14 

SCIM 

15 

SCIM 

16 

SCIM 

17 

SC/M 

18 

SCIM 

19 

SCIM 

20 

SCIM 

21 

SCIM 

22 

SCIM 

23 

SCIM 

24 

SCIM 

25 

SCIM 

26 

SCIM 

27 

SCIM 

28 

SCIM 

29 

SCIM 

30 

SCIM 

31 

SCIM 

32 

26), SCIM 

33 

SCIM 

34 

SCIM 

35 

SCIM 

36 

SCIM 

37 

SCIM 

38 

SCIM 

39 

SCIM 

40 

SCIM 

41 

SCIM 

42 

46 


c 


DECLARE  VARIABLES  FOR  BALLOON  DRIFT  TIHESTEP  SINULATION 
REAL  LATSTEP,LONGSTEP,DELTALAT,DELTALONC,NEULAT,NEWLONC 
INTEGER  TIME, STOPOAYS.TOTOAYS, PASS 
C  INITIALIZE  BALLOON  DRIFT  TIMESTEP  SINULATION  VARIABLES 
NEULAT  >  30.0 
PASS  •  0 
TOTOAYS  >  0 
STOPOAYS  •  362 

C  TINE  STEP  IS  1  HOUR:  24  TINE  INCREMENTS  PER  DAY 
STOPOAYS  <  (24  *  STOPOAYS) 


c 

FACTOR  FOR  RADIANS  TO  DEGREES 

SCIM 

43 

FAC  »  57.2957795 

SCIM 

44 

IUSYN>ICHAR('  ') 

SCIM 

45 

IF(NPOP.NE.O)  GO  TO  6 

SCIM 

46 

UPRE>0. 

SCIM 

47 

VPRE>0. 

SCIM 

48 

DUPRE<0. 

SCIM 

49 

OVPREsO. 

SCIM 

50 

6 

PQ*0. 

SCIM 

51 

DQ«0. 

SCIM 

52 

T0«0. 

SCIM 

53 

PRHsO. 

SCIM 

54 

DRH*0. 

SCIM 

55 

TRH<0. 

SCIM 

56 

URH*0. 

SCIM 

57 

VRH>0. 

SCIM 

58 

URH  s  0. 

SCIM 

59 

ua«0. 

SCIM 

60 

va«o. 

SCIM 

61 

paA>o. 

SCIM 

62 

OOA'O. 

SCIM 

63 

TOAiO. 

SCIM 

64 

UA>0. 

SCIM 

65 

VA«0. 

SCIM 

66 

PSH*0. 

SCIM 

67 

DSH<0. 

SCIM 

68 

TSN>0. 

SCIM 

69 

SPU  •  0. 

SCIM 

70 

SPV  »  0, 

SCIM 

71 

MONTH«MN 

SCIM 

72 

C 

PRESENT  latitude,  DEG 

SCIM 

73 

PHI  »  PHIR*FAC 

SCIM 

74 

C 

PRESENT  LONGITUDE,  DEG 

SCIM 

75 

THET  »  THETR*FAC 

SCIM 

76 

C 

PREVIOUS  LATITUDE,  DEG 

SCIM 

77 

PHM  *  PHI1R*FAC 

SCIM 

78 

c 

PREVIOUS  LONGITUDE,  DEG 

SCIM 

79 

c 

c 

196 

THET1  =  THET1R*FAC 

BEGINNING  OF  BALLOON  DRIFT  SIMULATION  DO  LOOP 

DO  LOOP  IS  FOR  THE  TOTAL  NUMBER  OF  DAYS  AS  INITIALIZED 

DO  197  K  »  1, STOPOAYS,! 

SCIM 

80 

C. . . . 

.FCORY  >  NORTH  COMPONENT  CORIOLIS  FACTOR  TINES  DISTANCE  FOR 

SCIM 

61 

c 

5  DEGREES  OF  LATITUDE 

SCIM 

82 

OY5  »  5000.*RI/FAC 

SCIM 

83 

DX5  •  DY5*COS(PHIR) 

SCIM 

84 

FCORY  =  DY5*SIN<PHIR)/(120.*FAC) 

SCIM 

85 

c.... 

IN  JACCHIA  OR  MIXED  ZONAL  MEAN-JACCHIA  HEIGHT  RANGE 

SCIM 

86 

8 

if(h.gt.90.0)  go  to  10 

SCIM 

87 

c.... 

.IN  A-0  DATA  HEIGHT  RANGE 

SCIM 

88 

IF  {H.LE.2S.0)  GO  TO  500 

SCIM 

89 

c.... 

IN  ZONAL  MEAN  OR  MIXED  ZONAL  MEAN  4D  HEIGHT  RANGE 

SCIM 

90 

GO  TO  200 

SCIM 

91 

c..,. 

IN  MIXED  JACCHIA'ZONAL  MEAN  RANGE,  NEED  TO  FAIR  DATA 

SCIM 

92 

10 

IF  (H.LT.IZO.)  GO  TO  20 

SCIM 

93 

c. . . . 

.FOLLOWING  IS  THE  PURE  JACCHIA  HEIGHT  RANGE  SECTION 

SCIM 

94 

c.... 

.JACCHIA  VALUES  AT  CURRENT  POSITION 

SCIM 

95 

CALL  JACCH(H,PHIR,THET,PH,DH,TN) 

SCIM 

96 

PHIN  >  PHIR  *  5.  /  FAC 

SCIM 

97 

THETE  «  THET  •  S. 

SCIM 

96 

c... 

..JACCHIA  VALUES  AT  CURRENT  POSITION'S  DEGREES  LAT,  FOR  DP/OY  AND 

SCIM 

99 

c 

OT/OV 

SCIM 

100 

CALL  JACCH(H.PHIN,THET,PHN.OHNJHN) 

SCIM 

101 

c,.. 

..JACCHIA  VALUES  AT  CURRENT  POSITION-S  DEGREES  LON,  FOR  DP/DX  AND 

SCIM 

102 

c 

DT/OX 

SCIM 

103 

CALL  JACCHCH, PHIR, THETE, PHE, ONE. THE) 

SCIM 

104 

c 

DP/DY  FOR  GEOSTROPHIC  WIND 

SCIM 

105 

0PY«PNN‘PH 

SCIM 

106 

c 

OP/DX  FOR  GEOSTROPHIC  WIND 

SCIM 

107 

DPXsPHE'PH 

SCIM 

108 

c 

OT/DX  FOR  THERMAL  WIND  SHEAR 

SCIM 

109 

DTX  «  THE  -  TH 

SCIM 

110 

c 

OT/DY  FOR  THERMAL  WIND  SHEAR 

SCIM 

111 

DTY  •  THN  •  TH 

SCIM 

112 

c 

CHANGE  NOTATION  FOR  OUTPUT 

SCIM 

113 

PGH»PH 

SCIM 

114 

DCHOH 

SOM 

115 

TGH«TH 

SCIM 

116 

CALL  WIND 

SCIM 

117 

UH  *  UGH 

SCIM 

116 

VH  a  VGH 

SCIM 

119 

H8  a  H  ♦  5. 

SCIM 

120 

CP  a  7.*PH/<2.*0H*TH) 

SCIM 

121 

CALL  JACCH(HB. PHIR, THEY, PB, OB. T8) 

SCIM 

122 

DTZ  a  (TB  «  TH)/S000. 

SCIM 

123 

c... 

..VERTICAL  MEAN  WIND 

SCIM 

124 

WGH  a  -CP*<UH*OTX/OX5  ♦  VH«OTY/DY5)/<G 

CP*OTZ  ♦  UN*OUH*VH*DVH) 

SCIM 

125 

c 

GO  TO  RANDOM  PERTURBATIONS  SECTION 

SCIM 

126 

GO  TO  600 

SCIM 

127 

c.. . 

.  FOLLOWING  IS  THE  MIXED  JACCHIA-ZONAL  MEAN  HEIGHT  RANGE  SECTION 

SCIM 

128 

c 

LOWER  HEIGHT  INDEX 

SCIM 

129 

20 

IHA  «  5«<INT(H)/S) 

SCIM 

130 

c 

UPPER  HEIGHT  INDEX 

SCIM 

131 

IH«  >  IHA  .  S 

SCIM 

132 

c 

lower  height  for  interpolation 

SCIM 

133 

HA  a  IHA*1. 

SCIM 

134 

c 

UPPER  HEIGHT  FOR  INTERPOLATION 

SCIM 

135 

HB  >  IH8*1. 

SCIM 

136 

c..., 

..JACCHIA  VALUES  AT  LOWER  HEIGHT,  CURRENT 

LAT -LON 

SCIM 

132 

CALL  JACCH<HA,PH1R, THET, PJA. OJA, TJA) 

SCIM 

138 

PHIN  a  PHIR  *  S.  /  FAC 

SCIM 

139 

THETE  a  THET  •  5. 

SCIM 

140 

c... 

..JACCHIA  VALUES  AT  LOWER  HEIGHT,  CURRENT 

LON-LAT.S  OECBEES 

SCIM 

141 

c 

LAT,  FOR  OP/OY  AND  OT/DY 

SCIM 

142 

CALL  JACCHCHA, PHIN, THET, PJN,DJN,TJN} 

SCIM 

143 

c... 

..JACCHIA  VALUES  AT  LOWER  HEIGHT,  CURRENT 

LAT-LON-5  DEGREES 

SCIM 

144 

c 

LON,  FOR  DP/DX,  AND  OT/OX 

SCIM 

145 

CALL  JACCH(HA, PHIR, THETE. PJE.OJE.TJE) 

SCIM 

146 

c 

JACCHIA  DP/OY  AT  LOWER  HEIGHT 

SCIM 

142 

OPXJAapJE-PJA 

SON 

146 

c 

JACCHIA  OP/DY  AT  LOWER  HEIGHT 

SCIM 

149 

OPYJAaPJN-PJA 

SCIM 

150 

c 

JACCHIA  OT/DX  AT  LOWER  HEIGHT 

SCIM 

151 

OTXJA  a  TJE  -  TJA 

SCIM 

152 

c 

JACCHIA  OT/DY  AT  LOWER  HEIGHT 

SCIM 

153 

OTYJA  a  TJN  •  TJA 

SCIM 

154 

c... 

..JACCHIA  VALUES  AT  UPPER  HEIGHT,  CURRENT 

LAT'LON 

SCIM 

155 

CALL  JACCH( HB, PHIR, THET, PJB.DJB.TJB) 

SCIM 

156 

PHIN  a  PHIR  ♦  5.  /  FAC 

SCIM 

157 

THETEaTHET-5 

SCIM 

158 

c... 

..JACCHIA  VALUES  AT  UPPER  HEIGHT,  CURRENT 

L0N-LAT*5  DEGREES 

SCIM 

159 

c 

LAT,  FOR  DP/OY  AND  OT/DY 

SCIM 

160 

CALL  JACCH(NB, PHIN, THET, PJN,OJN. TJN) 

SCIM 

161 

c... 

..JACCHIA  VALUES  AT  UPPER  HEIGHT,  CURRENT 

LAT-LON-S  DEGREES 

SCIM 

162 

c 

LON,  FOR  DP/OX  AND  OT/OX 

SCIM 

163 

CALL  JACCHCHB, PHIR, THETE, PJE.OJE.TJE) 

SCIM 

164 

c 

JACCHIA  OP/DX  FOR  GEOSTROPHIC  WINDS 

SCIM 

165 

48 


DPXJB  >  PJE  -  PJB  SCIN  166 

C  JACCHIA  OP/DY  FOR  CEOSTOPHIC  WINDS  SCIN  167 

DPYJB  >  PJN  ■  PJB  SCIN  168 

C  JACCHIA  DT/DX  FOR  THERMAL  WIND  SHEAR  SCIN  169 

OTXJB  >  TJE  -  TJB  SCIN  170 

C  JACCHIA  OT/Or  FOR  THERMAL  WIND  SHEAR  SCIM  171 

OTTJB  «  TJN  -  TJB  SCIN  172 

C _ ZONAL  MEAN  AT  LOWER  HEIGHT,  TO  BE  FAIREO  WITH  JACCHIA  SCIN  173 

CALL  GTERP(IHA,PHI,PCA,OGA,TGA,PG,OC,TG,OPrGA,OTTGA,UCA,UG)  SCIM  176 

C _ ZONAL  MEAN  AT  UPPER  HEIGHT,  TO  BE  FAIREO  WITH  JACCHIA  SCIM  175 

CALL  GTERP(IHB,PHI,PGB,OGB,TGB,PG,OG,TG,OPrGB,DTTGB,UGB,UG)  SCIM  176 

C . FAIRED  RESULTS  AT  LOWER  HEIGHT  SCIM  177 

IHSB  >  90  SCIM  178 

CALL  POTUVLPSP, DSP, TSP,USP,VSP, PHI, THET, IHSB, PSH, OSH, TSH,  SCIN  179 

t 


OPXSft.OPYSB.  OTXSB. 0TYS8«SPU,SPV) 

SCIM 

180 

PGA  *  PGA*(1.  ♦  PSH) 

SCIM 

181 

OGA  «  DGA*(1.  *  OSH) 

SCIM 

182 

TGA  ■  TGA*<1.  ♦  TSH) 

SCIM 

183 

PGB  «  PGB*(1.  ♦  PSH) 

SCIM 

184 

DGB  «  0GB*(1.  ♦  OSH) 

SCIM 

185 

TGB  ■  TG8*(1.  ♦  TSH) 

SCIM 

186 

UCA-UGA«SPU 

VCA<SPV 

UGB>UCB*SPU 

VGB>SPV 

DTXGA  >  OTXSB  *  TCA 

OTXGB  •  OTXSB  *  TGB 

OTYGA  •  TGA'DTYSB  ♦  DTYGA*(1.  ♦  TSH  *  OTYSB) 

OTYGB  =  TGB*DTYSB  ♦  DTYC8*(1.  ♦  TSH  ♦  DTYSB) 

OPXGA  >  DPXS8  •  PGA 

OPXGB  <  0PXS8  *  PG8 

DPYGA  >  PGA'OPYSB  *  0PYGA*(1.  *  PSH  «  OPYSB) 

OPYGB  >  PGB'OPYSB  *  OPYGB'd.  *  PSH  *  OPYSB) 

CALL  FAIRCPGA, OGA, TGA, PJA, OJA, TJA, IHA, PI, 01, T1, OPXGA, DPYGA, 
t  OPXJA, OPYJA,OPXA,OPYA,OTXGA, OTYGA, DTXJA,OTYJA,DTXA,OTYA) 

C . FAIREO  RESULTS  AT  UPPER  HEIGHT 

CALL  FAIR(PG8,0GB,TG8,PJB,0JB,TJB, IHB.PZ, 02, T2, OPXGB, OPYGB, 
i  0PXJ8,DPYJB,0PXB,0PYB, OTXGB, OTYGB, OTXJB, DTYJB,OTXB,OTYB) 

C . HEIGHT  INTERPOLATION  ON  FAIREO  P,0,T 

CALL  INTER2(P1,D1,T1,HA,P2,02,T2,H8,PH,0H,TH,H> 

C . HEIGHT  INTERPOLATION  ON  FAIREO  OP/OX, OP/OY 

CALL  INTERU(DPXA,DPYA,HA,OPXB,OPYB,HB,DPX,OPY,H) 

C . HEIGHT  INTERPOLATION  ON  FAIRED  OT/DX,DT/OT 

CALL  INTERW(DTXA,0TYA,HA,0TXB,0TTB,HB,DTX,0TY,H) 

C....  HEIGHT  INTERPOLATION  OF  WIND 

CALL  INTERU(UGA,VGA,HA,UG8,VGB,HB,USH,V$H,H> 
0USH:(UGB-UGA)/5000. 

DVSH=(VGB-VGA)/5000. 

C  CHANGE  OF  VARIABLES  FOR  OUTPUT 
PGH:PH 
OGH^OH 
TGH>TH 
CALL  WIND 
UHsUGH 
VH>VGH 

CP  «  7.*PH/(2.*0H*TH) 

DTZ  ■  (T2  -  T1)/5000. 

C . VERTICAL  MEAN  WIND 

WOH  •  -CP»(UH*0IX/0X5  ♦  VH*0TY/0Y5)/<G  *  CP»DT2  ♦  UH'DUH 
C  GO  TO  RANDOM  PERTURBATIONS  SECTION 

GO  TO  800 

C _  THE  FOLLOWIHG  SECTION  IS  FOR  ZONAL  MEAN  OR  MIXED  ZONAL  MEAN  40 

C  HEIGHTS 

C  UPPER  HEIGHT  INDEX 

200  IHG8  =  5*(1HT(H)/5)  ♦  5 

C  UPPER  HEIGHT 

HG8  •  IHG8*1. 

C _  ZONAL  MEAN  AT  UPPER  HEIGHT 
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CALL  GTERP(IHGB,PHI,PGB,OGB,TGB,PG,OG,TG,DPVGB,OTrGB,UGB.UG) 

IHSB  «  5«{IMT(M)/5)  ♦  5 
IF  (IHSB  .GT.  90>INSB  •  90 
C  UPPED  STATIONART  PERTURBATION  HEIGHT 

230  HSB  •  1HSB*1. 

C . STATIONARY  PERTURBATIONS  AT  UPPER  HEIGHT 

CALL  POTUV(PSP,OSP,TSP,U$P,VSP, PHI, THET, IHSB, PSB.DSB.TSB, 

$  DPXSB,DPYS8,  DTXSB,DTYSB,USB,VSB) 

C  LOUER  HEIGHT  INDEX 

IHGA  «  IHGB  -  5 
C  LOWER  HEIGHT  INDEX 

HGA  >  IHGA*1. 

C...  ZONAL  MEAN  AT  LOUER  HEIGHT 

CALL  GTERP(IHGA,PHI ,PGA,OGA,TGA,PG,DG,TG.OPYGA,OTYGA,UGA,UG) 
INSA<IHSB  -  5 

C  LOUER  STATIONARY  PERTURBATION  HEIGHT 

250  HSA  •  IHSA*1. 

C . STATIONARY  PERTURBATIONS  AT  LOWER  HEIGHT 

CALL  POTUV(PSP , DSP , TSP , USP , VSP , PH I , THE T , I HSA, PSA . OSA , TSA , 

S  OPXSA,OPYSA,  OTXSA,OTYSA,USA,VSA) 

CALL  INTERU(UGA,0.,HGA,UCB,0.,HGB,UGH,VGH,H) 

CALL  INTERU(USA,VSA, HSA, USB, VSB, HSB, SPU,SPV,H> 

USH^UGH-fSPU 

VSH>SPV 

0USH>((UGB-UGA)/(HGB-MCA)*(USB-USA)/(HSB-HSA))*.001 
OVSH> .  001*(VSB-VSA)/  ( HSB-  HSA ) 

C  FOR  MIXED  ZONAL  MEAN  -  40  SECTION 

IF(H.LT.30.0)  GO  TO  300 
C....  ZONAL  MEAN  VALUES  HEIGHT  INTERPOLATIONS 

CALL  INTER2(PGA,DGA,TGA,NGA,PGB,0GB,TGB,HG8,PGH,0GH,TCH,H> 

C . STATIONARY  PERTURBATION  HEIGHT  INTERPOLATION 

CALL  INTERZ<PSA,DSA,TSA,HSA,PSB,OSB,TSB,HSB,PSH,D$h,TSH,H> 

C  QUASI -BIENNIAL  VALUES 

CALL  QBOGEN 

C _  HEIGHT  INTERPOLATION  OF  ZONAL  MEAN  OP/DY  AND  OT/DY 

CALL  1NTERW(DPYGA,0TYGA,HGA,0PYCB,0TYGB,HCB,DPYG, 
t  0TYG,H) 

C . HEIGHT  INTERPOLATION  OF  STATIONARY  PERTURBATION  OP/DX  AND  OP/OY 

CALL  INTERW(DPXSA,DPYSA, HSA, 0PXSB,0PYSB, HSB, DPXS, OPTS, H> 

C . HEIGHT  INTERPOLATION  OF  STATIONARY  PERTURBATION  OT/DX  AND  DT/DT 

CALL  INTERW(DTXSA,DTYSA,HSA,DTXSB,0TYSB,H$B,DTXS,DTY5,H) 

C . UNPERTURBED  (MONTHLY  MEAN)  VALUES  FOR  OUTPUT 

TGH  »  TGH  •  (1.  ♦  TSH) 

PGH  »  PGH  •  (1.  ♦  PSH) 

DGH  •  DGH  *  (1.  *  OSH) 

C  TOTAL  OT/OX 

DTX  =  OTXS  *  TGH 

C  TOTAL  OT/DY 

OTY  »  TCH*OTYS  ♦  DTYG*(1.  ♦  TSH  ♦  DTYS) 

C  TOTAL  OP/OX 

OPX  =  DPXS  *  PGH 

C  TOTAL  DP/DY 

DPY  *  PGH»OPYS  ♦  0PYG*(1.  ♦  PSH  ♦  DPYS) 

C . UNPERTURBED  VALUES  PLUS  OBO  PERTURBATIONS 

PH  *  (1.  ♦  PO)  •  PGH 
DH  >  DGN  •  (1.  OQ) 

TH  =  (1.  ♦  TO)  •  TGH 
CALL  WIND 

C  GEOSTROPHIC  WIND  PLUS  QBO  WIND  PERTURBATIONS 

UH<UGH'HJO 
VH^VGHtVO 

CP  =  7.»PGH/(2.»0GH*TGH) 

DTZ  «  (TGB*{1.«TSB)  -  TGA*(1.*TSA))/5000. 

C . VERTICAL  MEAN  WIND 

C  THIS  CHANGE  WAS  MADE  5  AUG  92  DUE  TO  PHONE  CALL  WITH  MSFC  JEFFRIES 
C  UGH»-CP*(UGH*DTX/DX5*VGH*DTY/DY5)/(G*CP*DTZ«VGH*DUH»VGH*DVH) 

UGH«-CP*(UGH*DTX/DX5*VOH*OTY/OY5)/(G*CP*DTZ»UCH*OUM*VGH»DVH) 

C  GO  TO  RANDOM  PERTURBATIONS  SECTION 
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GO  TO  800 

C _  THE  FOLLOyiNG  1$  THE  MIXED  ZONAL  MEAN-40  SECTION 

C . GENERATE  GRID  OF  40  PROFILES  IF  PREVIOUS  HEIGHT  GE  30 

300  IF  (LOOK  .EQ.  DCALL  GEN40 
IHCK  •  24 
DO  310  KNO  •  1,3 
IKNO  >  IHCK  *  KNO 
IF  (IKN0.GT.26)IKND>26 
00  310  INO  •  1,4 
DO  310  JNO  •  1,4 

PCK(INO,JNO,KND)  <  P40(4*( INO- 1 )«JNO, IKNO) 

DCK(INO,JNO,KND)  «  D40(4*( INO- 1 >* JNO, IKNO) 

310  CONTINUE 

C . LAT-LON  INTERPOLATION  OF  40  DATA  AT  25  KM 

CALL  INTER4(  PHI,THET,2S,  P40,040,T4D,P4A,04A,T4A, 

S  DPX4,DPV4,DTX4,DTY4) 

C _  ZONAL  MEAN  PLUS  STATIONARY  PERTURBATIONS 

PB  >  PGB*(1.  *  PSB) 

C  P,0,T 

OB  >  DCB*(1.  *  DSB) 

TB  >  TGB*(1.  *  TSB) 

DPXB  <  PGB*DPXSB 

OPYB  •  PGB*OPYSB  *  DPYGB*(1.  *  PSB  *  OPYSB) 

DTXB  •  TGB*OTXSB 

DTYB  »  TGB*0TYS8  ♦  0TYGB*<1.  ♦  TS8  ♦  OTYSB) 

C _  HEIGHT  INTERPOLATION  BETWEEN  40  AT  25  AND  ZONAL  MEAN  AT  UPPER 

C  HEIGHT  OP/DX  AND  OP/OY 

CALL  I NTERU(0PX4 , DPY4 , 25 . , DPXB , OPYB , HSB , DPX , OPY , K ) 

C....  HEIGHT  INTERPOLATION  BETWEEN  40  AT  25  AND  ZONAL  MEAN  AT  UPPER 
C  HEIGHT  P,D,T 

CALL  INTER2(P4A,D4A,T4A,25.,PB,0B,TB,HG8,PGH,0GH,TCH,H) 

C _  HEIGHT  INTERPOLATION  BETWEEN  40  AT  25  AND  ZONAL  MEAN  AT  UPPER 

C  HEIGHT  DT/OX  AND  DT/OY 

CALL  INTERW(DTX4,DTY4,25.,0TXB,DTY8,HSB,DTX,0TY,H) 

IF  (IOPO.Ea.2)  GO  TO  350 
C  QUASI  BIENNIAL  PERTURBATIONS 

ULL  QBOGEN 

C  ADO  QBO  PERTURBATIONS  TO  P,D,T 

350  PH»PGH*(1.*P0) 

OH>OGN*(1.*OQ) 

TH»TGH*<1.*TQ) 

CALL  WINO 

C  ADD  QBO  WINO  PERTURBATIONS 

UHiUGH^UO 
VH»VGH*VO 

CP  «  7.*PGH/(2.*0GH»TGH) 

DTZ  »  (TB  -  T4A)/(1000.*(HG8  -  25.)) 

C . VERTICAL  MEAN  WIND 

WCH»-CP*(UGH»DTX/DX5*VGH*OTY/OY5)/(G*CP*DTZ*UGM*DUH*VGM*OVH) 

C  GO  TO  RANDOM  PERTURBATIONS  SECTION 

2000  FORNAT('  4-0  DATA  AFTER  ADJUSTMENTS'/'  LATITU0E'/3X, 16F8.3) 

2001  FORMAT('  LONCITUDE'/3X,16F8.3) 

2007  FORMAT!'  PRESSURE') 

2002  FORMAT(1X,I2,16F8.0) 

2003  FORMAT!'  DENSITY') 

2005  FORMAT!'  TEMPERATURE') 

2004  F0RMAT!1X,I2,16F8.5) 

2006  FORMAT!1X,I2,16F8.2) 

GO  TO  800 

500  IF  (H.GE.0.0)  GO  TO  510 
IF  (H.LT. -0.015)  GO  TO  505 
C  IF  -15  METER  LE  N  LT  0  ,  H  IS  SET  TO  0 
H  •  0. 

GO  TO  510 

C  NO  MORE  COMPUTATIONS  TO  BE  MADE  IF  HEIGHT  LT  -5  M 
505  NMORE  >  0 
RETURN 

C . GENERATE  GRID  OF  40  PROFILES  IF  PREVIOUS  HEIGHT  GE  30 
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510  IF  (LOOK  .EQ.  DCALL  GEN40 
C  LOWER  HEIGHT  INDEX 

IHA<INT(H) 

C  LOWER  HEIGHT  INDEX 

HA  •  IHA*1. 

IWSX  >  IWSYM 
IHCK>1HA-1 
DO  511  KNO*1,3 
UNO  >  IHCK  4^  KND 
IF  (IKNO.LT.DUND  >  1 
IF  (IKNO.GT.26)UNO  >  26 
DO  511  IND>1,4 
DO  511  JNO  •  1,4 

PCK(IN0,JN0,KN0}*P40(4*(IN0-1)4'JND,UN0) 

DCK(IN0,JN0,KN0)440(4*(IN0-1)*JND,UND) 

511  CONTINUE 
C  UPPER  HEIGHT  INDEX 
IHB  «  IHA  *  1 
IF(INB.LE.25>  GO  TO  513 
IHA*24 
HA>24. 

IHB°25 

C  UPPER  HEIGHT 

513  HB  «  IHB*1. 

C . LAT-LON  INTERPOLATION  OF  40  VALUES  AT  UPPER  HEIGHT 

515  CALL  INTER4(  PHI ,THET, IHB,  P4D,D40,T40,PB,DB,TB, 

S  0PX4B,0PY4B,DTX4B,0Ty4B) 

IF(IHA.Ea.O.ANO.PB*OB*TB.LE.O.)GO  TO  520 
GO  TO  540 
520  IHB<IHB*1 

C . LOOP  TO  FIND  LOWEST  VALID  HEIGHT 

HBiHB-^l. 

GO  TO  515 

540  IF(IHA.GT.0)CALL  INTER4(  PHI ,THET, IHA,  P4D,D4D,T4D, 

$  PA ,DA , TA , 0PX4A , OPYAA , DTX4A , 0TY4A ) 

IFdWSYM.EQ.ICHARC*'))  IWSX<IUSYH 

IF(IHA.Ea.0.OR.(PA*DA*TA.LE.0.AND.IHA.LT.10.AN0.P8*D8*T8.GT.0.)) 
1GO  TO  550 
GO  TO  600 

C . LAT-LON  INTERPOLATION  OF  40  VALUES  AT  LOWER  HEIGHT 

550  CALL  INTER4(  PHI,THET,0,  P40,D4D,T4D, 

$  PA,0A,TA,0PX4A,DPY4A,DTX4A,0TY4A) 

IF(  IWSYN.EQ.  ICHARC  '*'  ))  IUSX=IWSYN 
IF(TA-TB)560,570,560 
560  IF(TA»TB.LE.0.0)  GO  TO  570 
TZ  •  (TA-TB)  /  ALOC(TA/TB) 

GO  TO  575 
570  TZ«TA 

C  ...COMPUTES  HEIGHT  OF  SURFACE 

575  HA  <  HB 
IF(PB*PA.LE.O.O)GO  TO  576 

HA  >  HB  *  0.28705*TZ*ALOG(PB/PA)/C 

576  IF(H.GT.HA  -  0.04)GO  TO  600 
PH:0. 

DH>0. 

TH»0. 

PGH'O. 

DGH<0. 

TGH«0. 

GO  TO  800 

C . HEIGHT  INTERPOLATION  OF  P,0,T 

600  CALL  INTER2(PA,DA,TA,HA,PB,DB,TB,HB,PGH,DGH,TGH,N) 

C . HEIGHT  INTERPOLATION  OF  OP/DX  AND  DP/DY 

CALL  INTERU(DPX4A,DPY4A,HA,0PX48,DPY4B,HB,0PX,DPY,H) 

C . HEIGHT  INTERPOLATION  OF  DT/DX  AND  DT/DY 

CALL  INTERU(OTX4A,DTY4A,HA,OTX4B,OTY4B,NB,OTX,DTY,N) 

C  CHANGE  OF  NOTATION  FOR  OUTPUT 
PH  •  PGH 
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DN  >  DGN 
TN  •  TCM 

IF<PH*OM*TM.LE.O.)  GO  TO  800 
CALL  WIND 

C  CHANGE  OF  NOTATION  FOR  OUTPUT 

UH  >  UGH 
VH  >  VGH 

CP  •  /.♦PCH/(2.*DCH*TCH) 

DTZ  >  (TB  -  TA)/(1000.*(HB  -  HA)) 

C....  VERTICAL  MEAN  UINO 

UGH  *  -CP*(UGH*DTX/DX5  ♦  VGH*OTY/Or5)/(G*CP*DTZ*UH*OUH*VH*OVH) 
C  OBO>0  IF  H  LT  10 

IF  (H.LT.10.)  GO  TO  800 
IF  (lOPQ.EQ.Z)  GO  TO  6S0 
C  COMPUTES  GUASI  BIENNIAL  PERTURBATIONS 
CALL  QBOGEN 

C  ADDS  OBO  PERTURBATIONS  TO  P,0,T 
650  PH*PGH*(1.*PO) 

DH>0GH*(1.*00) 

TH«TGH*(1.*TO) 

C  ADDS  OBO  UINO  PERTURBATIONS  TO  U,V 
UH^UGH^UO 
VH«VGH*Va 

C . THE  FOLLOUING  IS  THE  RANDOM  PERTURBATIONS  SECTION 

C . NO  RANDOM  PERTURBATIONS  IF  lOPR  GT  1 

800  CONTINUE 

IF(H  .GE.  30.)GOTO  512 
IF(IPRT.NE.O)GOTO  512 
URITE(6,2000)  (GLAT( I ), I<1 ,NG) 

URITE(6,2001)  (GL0N(I),I°1,NG) 

URITE(6,2007) 

DO  504  I<1,26 
IH«I-1 

URITE(6,2004)IH,(SP4(J,I),J*1,NG) 

URITE(6,2002)  lH,LP«)LJ,n,)*1,HG) 

504  CONTINUE 

URITE(6,2003) 

DO  507  I  •  1,26 
IH  •  I  -  1 

URITE(6,2004)IH,(S04(J,I),J°1,NG) 

507  URITE(6,2004)IH,(D40(J,I),J>1,NG) 

URITE(6,2005) 

DO  506  I  >  1,26 
IH  •  I  -  1 

URITE(6,2004)IN,(ST4(J,I),J=1,NG) 

506  URITE(6,2006)IH,(T40(J,I),J=1,HG) 

IPRT«IPRT*1 
512  CONTINUE 

IF(NPOP.EO.O)GO  TO  840 
IF  (IOPR.GT.1)  GO  TO  830 

C . INTERPOLATES  RANDOM  UINO  MAGNITUDES  TO  HEIGHT  H,  LATITUDE  PHI 

CALL  INTRUV(UR,VR,N,PHI,SUN,SVH) 

CALL  INTR25(PLP,DLP,H,PHI,PLPH,DLPH) 

CALL  INTR25(TLP,DLP,H,PHI,TLPH,DLPM) 

CALL  I NTR25(ULP , VLP , H , PH  I ,ULPH , VLPH ) 

CALL  INTR25(UDL,VDL,H,PHI,UOL2,VDL2) 

CALL  INTR25(UDS,VDS,H,PNI,UDS2,VDS2) 

CALL  INTRU(UR,H,SUH) 

SUHL«SORT(ULPH*ABS(SUH}) 

SUHS'SORT ( ( 1 . -ULPH )*ABS ( SUH ) ) 

SVHL>SORT(VLPH*ABS(SVH}) 

SVHS*SORT ( ( 1 . ■ VLPH )*ABS( SVH ) ) 

SUH  >  SORT(ABS(SUH)) 

SVH  •  SORT(ABS(SVH)) 

IF(H.CE.25.)G0T0  805 

C .  IF  H  LE  20  USE  40  DATA  RANDOM  P,D,T  SIGMAS 

IF(N.LE.20.)GOTO  810 

C . INTERPOLATE  PR,0R,TR  ARRAYS  TO  GET  P,D,T  SIGMAS  AT  HEIGHT  H, 


SCIM  436 
SCIM  437 
SCIM  438 
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c 


LATITUDE  PHI 

CALL  IITERP(2$.,PHI,PR,DR,TR,SPHG,S0NG,STHC) 
GO  TO  810 
80S  CONTINUE 


CALL  RTERP(H,PHI,PR,DR,TR,SPH,SON,STH) 

GO  TO  820 

C . LAT-LON  INTERPOLATION  ON  P,0,T  SIGMAS  AT  LOWER  HEIGHT 

810  CALL  INTERAC  PHI.THET, IHA,  SPA,S04,ST4,PA,DA.TA, 

$  DPX.DPY.DTX.DTY) 

C .  LAT-LON  INTERPOLATION  ON  P,D,T  SIGMAS  AT  UPPER  HEIGHT 

CALL  INTER4(  PHI.THET, IH6,  SP4.S04,ST4,PB,0B,TB, 

%  DPX.DPY.DTX.DTY) 

C . HEIGHT  INTERPOLATION  OF  SIGMAS 


CALL  INTERZ(PA,OA,TA,  HA.PB.OB.TB,  HB.SPH.SOH.STH.H) 
IF(PH.LE.O.O.OR.DH.LE.O.O.OR.TH.LE.O.O)GO  TO  82S 
IF(H.LE.20.)GOTO  820 
FH  *  1.  -  0.2»(25.  -  H) 

SPH  »  FH*SPHG  +  <1.  -  FH)*SPH 
SON  »  FH*SOHG  +  (1.  -  FH)*SOH 
STH  »  FH*STHG  ♦  <1.  -  FH)*STH 

C . HEIGHT  DISPLACEMENT  BETWEEN  PREVIOUS  AND  CURRENT  POSITION 

820  DZ  >  HI  -  H 

SPHL:SaRT(PLPN*ABS(SPH)) 

SPHS>SORT( ( 1 . -PLPH )*ABS(SPH ) ) 

SOHL>SORT(OLPH*ABS(SOH)) 

SOHS:SORT( ( 1 . -DLPH )*ABS(SDH ) ) 

STHL>SaRT(TLPH*ABS(STH)) 

STHS>SORT( ( 1 . -TLPH)*ABS(STH ) ) 

SPH  >  SORT(ABS(SPH)) 

SDH  >  SORT(ABS(SOH)) 

STH  >  SORT(ABS(STH)) 


C . COMPUTES  HORIZONTAL  DISPLACEMENT  DX  BETWEEN  PREVIOUS  AND  CURRENT 

C  POSITION,  HORIZONTAL  SCALE  HL,  AND  VERTICAL  SCALE  VL 

C _  COMPUTES  PERTURBATION  VALUES  PRH,DRH,TRH,URH,VRH  AND  WRH 

CALL  PERTRB 

C  ADDS  RANDOM  PERTURBATIONS  TO  PH.DH.TH 

PH  «  PH*<1.  ♦  PRH) 

OH  «  OH*(1.  ♦  DRH) 

TH  >  TH*<1,  ♦  TRH) 

C  ADOS  RANDOM  WINDS  TO  UK,VH,WN 

UHMJHtURN 
VH>VH*VRH 
WH^WGH^WRH 

C . SETS  PREVIOUS  RANDOM  PERTURBATION  IN  P,D,T  TO  CURRENT 

C  PERTURBATIONS,  FOR  NEXT  CYCLE 


825  RP1S>  PRHS 
R01S>  ORHS 
RT1S:  TRHS 
RP1L»PRHL 
RDU^RHL 
RT1L=TRHL 

C _  SETS  PREVIOUS  MAGNITUDES  TO  CURRENT  VALUES,  FOR  NEXT  CYCLE 

SP1S>SPHS 

S01S<  SONS 

ST1S-STHS 

SPIL^SPHL 

S01L>S0HL 

STIL-STHL 


C . SETS  PREVIOUS  WIND  PERTURBATION  VALUES  TO  CURRENT  VALUES, 

C  FOR  NEXT  CYCLE 


RUIS^HS 

RV1S>VRHS 

RU1L*URHL 

RV1L«VRHL 

RW1<WRH 


C . SETS  PREVIOUS  WIND  PERTURBATION  MAGNITUDES  TO  CURRENT  VALUES, 

C  FOR  NEXT  CYCLE 

SU1S°SUHS 


SCIM  504 
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SCIM  506 
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SCIM  512 
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SCIM  515 
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SCIM  517 
SCIM  518 
SCIM  519 
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SCIM  521 
SCIM  522 
SCIM  523 
SCIM  524 
SCIM  525 
SCIM  526 
SCIM  527 
SCIM  528 
SCIM  529 
SCIM  530 
SCIM  531 
SCIM  532 
SCIM  533 
SCIM  534 
SCIM  535 
SCIM  536 
SCIM  537 
SCIM  538 
SCIM  539 
SCIM  540 
SCIM  541 
SCIM  542 
SCIM  543 
SCIM  544 
SCIM  545 
SCIM  546 
SCIM  547 
SCIM  548 
SCIM  549 
SCIM  550 
SCIM  551 
SCIM  552 
SCIM  553 
SCIM  554 
SCIM  555 
SCIM  556 
SCIM  557 
SCIM  558 
SCIM  559 
SCIM  560 
SCIM  561 
SCIM  562 
SCIM  563 
SCIM  564 
SCIM  565 
SCIM  566 
SCIM  567 
SCIM  568 
SCIM  569 
SCIM  570 
SCIM  571 
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SV1S<SVHS 

SUU>SUML 

SV1L>SVNL 

SU1*SUH 

C . SETS  PREVIOUS  HEIGHT  TO  CURREHT  HEIGHT,  FOR  NEXT  CYCLE 

830  H1  <  H 

C . SETS  PREVIOUS  LATITUDE  TO  CURRENT  LATITUDE,  FOR  NEXT  CYCLE 

PNI1R>PHIR 

C . SETS  PREVIOUS  LONGITUDE  TO  CURRENT  LONGITUDE,  FOR  NEXT  CYCLE 

THET1R>THETR 

C  SETS  NFtORE  TO  COMPUTE  MORE  DATA  ON  NEXT  CYCLE 
840  NMORE  <  1 

C . NO  MORE  DATA  IF  P,  0,  OR  T  LEO  0 

IF(PH»DH*TH.LE.O.)  RETURN 
CALL  STDATM(H,TS,PS,DS) 

IF  ((PS*OS*TS).GT.O.)  GO  TO  870 
PGHP>0. 

DGHP>0. 

TGHP«0. 

PHP»0. 

DHP>0. 

THP«0. 

GO  TO  880 

870  PGHP*100.»(PGH-PS)/PS 
OGHP«100.*(OGH-OS)/OS 
TGHP*100.*(TGH-TS)/TS 
PHP«100.*(PH-PS)/PS 
0HP*100.*(0H-DS)/0S 
THP*100.*(TH-TS)/TS 
C  CONVERTS  080  P,0,T  TO  PERCENT 
880  PO*100.*PO 
OQ>100.*00 
T0»100.*T0 

C  CONVERTS  RANDOM  P,D,T  TO  PERCENT 
PRH*100.*PRH 
0RH>100.*DRH 
TRH<100.*TRH 
PRHS<iaO.*PRHS 
DRHS<100.*ORNS 
TRHS<100.*TRHS 
PRHL*100.*PRHL 
0RHL<100.*DRHL 
TRHL>100.*TRHL 
SPHS  <  100.*SPHS 
SONS  >  100.*SONS 
STHS  •  100.»STHS 
SPHL  s  TOO.«SPHl 
SOHL  •  100.*SOHL 
STHL  «  lOO.'STHL 

C  CONVERTS  WIND  SHEAR  TO  M/S/KH 
DUH  >  DUN  ■  1000. 

OVH  >  DVH  •  1000. 

PQA»PQA»100. 

DQA^QA'IOO. 

TQA=TQA»100. 

SPH<SPH*100. 

SOH*SOH*10O. 

STH»STH*100. 

PSH>PSH*100. 

DSH>OSH*iaO. 

TSH>TSH*100. 

IF(NPOP.NE.O)  GO  TO  920 

UPRE*UGH 

VPRE'VGH 

0UPRE4UH/1000. 

DVPRE*OVH/1000. 

RETURN 

920  IF  (lOPP.NE.O) 
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*  URITE(10PP,9S1)H,PHI,TNET,DGHP,TGH,UGH,VGH,WGH,SDHL,STNL,  SCIM  640 

t  SUHL.SVHL  SCIN  641 

951  FORNAT(F5.1,7F7.2,4F6.2)  SCIM  642 

C  URITE(6,900)  H,PH1 ,THET,PGH,DGH,TGN,UGH,CHAR(IWSYM),  SCIN  643 

C  1  VGH,PH,DH,TH,UH,CHAR(1USYM),VH,DUH,  SCIN  644 

C  S  OVH,$UH,IET,PCHP,DGHP,rGHP,UGH,PHP,DHP,THP.WH.PSH,DSH,TSH,  SCIN  645 

C  S  SPU,SPV,PQ,DO,Ta,UO,  SCIM  646 

C  S  VO,PQA,DaA,TaA,UA,VA,PRHS,DRHS,TRHS,URHS,VRHS,SPHS,SOHS,STHS.  SCIM  647 

C  1$UHS.SVHS,PRHL,ORNL,TRHL,URHL, VRML, SPNL.SOHL.STHL, SUHL.SVHL,  SCIN  648 

C  2PRH,DRH,TRH,URH,VRH,SPH,SOH,STH,SUH,SVH  SCIM  649 

C900  FORNAT(1X,F6.2,2F7.2,2(2E9.3,2F6.0,A1.F5.0),2F5.1,23X,F6.2/1X,  SCIM  650 

C  1  I5,14X,2(F8.1,'  '),F6.1,'  ',E10.2,1X,  SCIM  651 

C  &  2<F8.1,'  •),F6.1,'  ',F10.2,11X,  SCIM  652 

C  23F5.1,2F5.0,'  SP'/102X,3F5. 1 ,2F5.0, '  O8O'/102X,3F5.1 ,2F5.0, '  MAG'/SCIM  653 

C  3  102X,3F5.1,2F5.0.'  RAHS'/102X,3F5.1.2F5.0,'  SIGS'/  SCIM  654 

C  4102X,3F5.1,2F5.0,'  RANL',/  SCIM  655 

C  5102X,3F5.1,2F5.0,'  SIGL',/  SCIM  656 

C  6102X,3F5.1.2F5.0.'  RAHT',/  SCIM  657 

C  7102X,3F5.1,2F5.0.'  SIGT',/)  SCIM  658 


C  NEXT  THREE  LINES  PRINT  OUT  LAT.LONG,  E-W  WIND,  N-S  UINO 
C  PRINT*, 'LAT,  LONG,  E-W  UINO,  N-S  WIND’ 

C  URITE(6,9091)  PHI ,THET,UH,CHAR( IWSYM),VH 
C9091  FORMATE  IX, 2(F7. 2, 6X),F6.0,A1,F5.0,/) 

C  MULTIPLY  WIND  SPEED  AND  WIND  DIRECTION  BY  TIME  STEP  TO 
C  CREATE  WIND  VECTORS  FOR  THE  TIME  STEP  IN  KILOMETERS 

C  GIVEN  A  TIME  STEP  OF  1  HOUR  WITH  60  SEC/MIN,  60  MIN/HR 

LONGSTEP  »  UH  •  60  •  60  •  1  /  1000 
LATSTEP  »  VH  •  60  *  60  *  1  /  1000 

DELTALAT  =  LATSTEP  /111 

DELTALONG  >  LONGSTEP  /  ( (COSEPHI /FAC) )  *  111.4) 

OELTALONG  «  MQOEOELTALONG, 360.0) 

C  APPLY  WIND  VECTORS  TO  LAT/LONG  (COMPUTE  NEW  LAT/LONG) 

IF  (ABS(NEWLAT).EO.  90)  THEN 
NEWLAT  :  PHI  -  DELTALAT 
ELSE 

NEWLAT  :  PHI  *  DELTALAT 
END  IF 

IF  (NEWLAT  .GT.  90)  THEN 
NEWLAT  <  90-(NEWLAT  -  90) 

ELSE 

CONTINUE 
END  IF 

IF  (NEWLAT  .LT.  -90)  THEN 
NEWLAT  <  ABS(NEWLAT)  -  180 
ELSE 

CONTINUE 

ENDIF 

PHI1R  «  PHI  *  FAC 
PHIR  «  NEWLAT  *  FAC 
PHI  <  NEWLAT 
NEWLONG  <  THET  -  OELTALONG 
IF  (NEWLONG  .GT.  180)  THEN 
NEWLONG  <  NEWLONG  -  360 
ELSE 

CONTINUE 

ENDIF 

IF  (NEWLONG  .LT.  -180)  THEN 
NEWLONG  :  NEWLONG  *  360 
ELSE 
CONTINUE 
ENDIF 

THET1R  •  THET  •  FAC 
THETR  <  NEWLONG  *  FAC 
THET  •  NEWLONG 

C  URITE(6,9092)  LATSTEP, LONGSTEP, DELTALAT, OELTALONG, 

C  SNEWLAT, NEWLONG 

C9092  FORMAT(1X,'LATSTEP(KM)  LONGSTEP(KM)' ,1X, 

C  S'DELTALAT(DEG)  DELTALONG(OEG)  NEWLAT(OEG)  NEWLONG(DEG)',/, 
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C  *1X,f5.0,8X,f6.0,7X,F5.0,9X,F6.0,12X,F5.0,7)(,F6.0,/) 

C  CHECK  LONGITUDE  PICKET;  RECORD  LATITUDE  IF  LONGITUDE  IS  MET 

C  CURRENT  WRITE  STATEMENT  WRITES  LATITUDE  TO  A  FILE 

IF  (NEWLONC.lt. -42.0  .AND.  NEWLONG.GT . -46.0  .ANO.  PASS.EO.O) 
PASS  >  1 

C  WRITE(9,9093)  NEWLAT 

C9093  FORNAT(1X,F6.2) 

ELSE 

CONTINUE 

ENDIF 

C  INCREMENT  TIME  ANO  CHECK  STOPPING  RULE 

TIME  •  TIME  ♦  1 
IF  (TIME  .GT.  23)  THEN 
DATCOUNT  <  DAYCOUNT  *  1 
C  PRINT*, 'INSIDE  TIME  LOOP' 

TIME  •  0 

C  WRITE(6,9094)  LATSTEP.LONGSTEP.DELTALAT.OELTALONG, 

C  SNEWLAT.NEWLONG 

C9094  FORMAT(1X,'LATSTEP(KH)  LONGSTEP(KM)' , IX, 

C  S'DELTALAT(DEG)  OELTALONG(OEG)  NEWLAT(DEG)  NEWLONG(DEC)' ,/, 

C  $1X,F5.0,8X,F6.0,7X,F5.0,9X,F6.0,12X,F5.0,7X,F6.0,/) 

IF  (DAYCOUNT  .GT.  30)  THEN 
C  PRINT*, 'INSIDE  DAYCOUNT  LOOP' 

MN  <  MN  *  1 
MONTH  =  MN 

C  PRINT*, 'MONTH  IS  CHANGED' 

DAYCOUNT  :  0 
IF  (MONTH  .GT.  12)  THEN 
lYR  =  lYR  ♦  1 
MN  >  1 
MONTH  :  MN 

C  WRITE(6,909S)  lYR 

C  9095  FORHAT(1X,I2,/) 

ELSE 

CONTINUE 

ENDIF 

C  WRITE(6,9099)  MONTH 

C  9099  F0RMAT(1X,I3,/) 

ELSE 

CONTINUE 

ENDIF 

ELSE 

CONTINUE 

ENDIF 

C  LOCATION  AFTER  10  DAYS  OF  24  HOURS  PER  DAY  TIME  INCREMENTS 

IF  (TOTDAYS  .ED.  240)  THEN 
WRITE (7, 9096)  NEWLAT, NEWLONG 

9096  FORMAT(1X,F6.2,2X,F7.2) 

ELSE 

CONTINUE 

ENDIF 

C  LOCATION  AFTER  30  DAYS  OF  24  TIME  INCREMENTS  PER  DAY 

IF  (TOTDAYS  .EQ.  720)  THEN 
WRITE(9,9097)  NEWLAT, NEWLONG 

9097  FORMAT(1X,F6.2,2X,F7.2) 

ELSE 

CONTINUE 

ENDIF 

C  LOCATION  AFTER  360  DAYS  OF  24  TIME  INCREMENTS  PER  DAY 

IF  (TOTDAYS  .ED.  8640)  THEN 
URITE(8,909a)  NEWLAT, NEWLONG 

9098  FORMAT(1X,F6.2,2X,F7.2) 

ELSE 

CONTINUE 
ENDIF 

TOTDAYS  •  TOTDAYS  ♦  1 

IF  (NEWLONG  .LT.  -46.0  .AND.  NEWLONG  .GT.  -50.0)  THEN 
PASS  •  0 


THEN 


C 

C 


C  ELSEIF  (NEULONG  .LT.  -38.0  .AND.  NEULONG  .GT.  -42.0)  THEN 

C  PASS  •  0 

C  ELSE 

C  CONTINUE 

C  ENOIF 

197  CONTINUE 

C  197  PRINT*, 'END  OF  DO  LOOP' 

RETURN 

END 


SC  IN  639 
SCIM  660 
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